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Abstract 

A parton-shower model for hadronic collisions based on the emission properties of QCD 
dipoles is proposed. This proposal therefore extends the well-known radiation pattern of 
pure final-state colour dipoles to QCD initial-state radiation, both of which are treated 
perturbatively. Corresponding dipole splitting functions are derived and the kinematics 
of all dipole splittings is discussed. Application to hadron production in electron-positron 
annihilation, to Drell-Yan lepton-pair and QCD jet production yields encouraging results. 
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1 Introduction 

In the past decades, Monte Carlo event generators including QCD parton-shower routines, such 
as PYTHIA [1—5], HERWIG [6,7], or ARIADNE [8] have been very successful in correctly describ- 
ing, both qualitatively and quantitatively, a large range of QCD-related phenomena at different 
colliders, at different energies, and with different initial states. The success of these programs is 
based on good approximations in their treatment of logarithmically enhanced emission of QCD 
particles in soft and/or collinear regions of phase space. 

In conventional parton showers such as the ones in PYTHIA [9-11] and HERWIG [12], this is 
achieved by an expansion around the collinear limit. This manifests itself in the ordering of 
subsequent emissions by virtual masses supplemented with an explicit veto on increasing emis- 
sion angles or by an ordering by emission angles, respectivelyEl Alternatively, perturbative QCD 
cascades can be formulated in terms of splitting colour dipoles rather than partons. This has 
been realized in the shower algorithm in ARIADNE [8], which is based on the Colour Dipole 
Model (CDM) [15-17]. Splitting the dipoles and ordering the emissions in relative transverse 
momenta of subsequent splittings is equivalent to an expansion around the soft limits of the 
radiation process. In [15] it has been argued that such a dipole shower quite naturally fulfils 
the requirements of quantum coherence, which, for the parton showers, lead to angular ordering 
of subsequent emission, see e.g. [18]. It is interesting to note that the ARIADNE shower yields 
results, which show a similar or even better agreement with data from electron-positron anni- 
hilation into hadrons [19-23]. However, in the CDM initial-state radiation (ISR), i.e. parton 
emission off incoming partons, is not treated explicitly but taken into account by redefining ISR 
as final-state radiation (FSR) off hadron remnants [24] . To correctly model ISR in this picture, 
non-perturbative corrections have to be applied, cf. sec. (2) Equipped with such non-perturbative 
components in its modelling of initial-state associated radiation, ARIADNE also succeeded in de- 
scribing a wealth of DIS data in a very reassuring way, see for instance [25]. To some extent, the 
reason for this excellent performance in describing e + e~ and DIS data is not entirely understood. 
The cause could be better treatment of small- x effects in the DIS case, which are assumed to be 
an issue also for the forthcoming LHC. Equally well, it could be just the effect of a careful tuning 
of the additional non-perturbative parameters in the case of DIS. Another idea is related to a 
supposedly improved simulation of single, potentially non-global, potentially large logarithms 
stemming from soft corners of emission phase space. This appears as a consequence of the fact 
that the leading 1/Nq terms of such contributions are better accounted for if the description of 
the radiation is based on the dipole structure in both matrix element and phase space [26] . This 
blurred picture of, on the one hand, delivering overwhelming agreement with data of various 
measurements and, on the other hand, lacking clear determination of the reason for this success 
provides a fair, but not the only motivation for trying out an alternative path in modelling ISR 
arising from colour dipoles. 

In view of the upcoming LHC era, Monte Carlo event generators are undergoing an intensive 
overhaul, leading essentially to complete rewrites of the codes [27-30] or to the construction of 
entirely new programs such as SHERPA [13] in the modern, object-oriented programming lan- 
guage C++. Apart from issues related to maintenance, a number of improvements concerning 
physics simulation motivated the construction of new event generators. First of all, the shower 
algorithms themselves, forming an essential part of the event generators, have been improved: 
in PYTHIA, a k± ordered parton shower has been introduced [4,31] in order to better account for 
coherence effects. There is also a dramatically extended model of multiple parton interactions. 



The current parton shower implementation of SHERPA [13], APACIC++ [14], is very similar to the well- 
established virtuality-ordered PYTHIA shower. 
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In HERWIG, a new formulation of angular ordering [32] better embeds Lorentz invariance and 
provides an improved treatment of those regions, where the original HERWIG shower over- or 
undercounted parton emissions. In addition, a new parton-shower formulation has been devel- 
oped based on Catani-Seymour dipole factorization [33-36], and steps have been undertaken 
in the development of yet another QCD shower formulation, which uses antenna functions [37]. 
For all these recent developments, a common denominator has been to put more emphasis on 
the notion of a colour-connected partner of the splitting parton and thus a reduction of the 
difference between parton and dipole showers. Especially, the showers based on either Catani- 
Seymour [38,39] or antenna subtraction kernels [40,41], aim at an improved matching with exact 
higher-order QCD matrix elements. In fact, considering the need for increased precision, this 
systematic inclusion of higher orders in the perturbative expansion of QCD has been a dramatic 
and recent improvement of the paradigm underlying building and using multipurpose Monte 
Carlo event generators. 

In the matching approach, the exact next-to-leading order matrix-element result is consistently 
combined with the resummation of the parton shower [42-45] such that the overall result cor- 
rectly reproduces the corresponding NLO total cross section and the first additional hard QCD 
emission. This has been first implemented for specific processes in MC@NL0 [46] on the basis of the 
Frixione-Kunszt-Signer subtraction [47]. This method depends to some extent on the details of 
the parton shower and also has some residual dependence on the process in question. With the 
PDWHEG-approach a shower-independent matching solution [44, 45] extends the original MCSNLO 
proposal and the appearance of negative weights, which are present in the former method, can 
be circumvented. 

In alternative approaches, sequences of tree- level multileg matrix elements with increasing final- 
state multiplicity are merged with the parton shower to yield fully inclusive samples correct 
at leading logarithmic accuracy by avoiding double-counting and missing phase-space regions. 
A first approach, known as the CKKW merging approach, has been presented for the case 
of electron-positron annihilations into jets in [48]; later it has been extended to hadronic colli- 
sions [49] and it has been reformulated to a merging procedure - called LCKKW - in conjunction 
with a dipole shower [50]. A further method, the MLM method, has been developed also aiming 
at a merging of matrix elements and parton showers. It uses a different way in generating the 
inclusive samples based on a geometric interpretation of the full radiation pattern in terms of 
cone jets [51,52]. These different algorithms have been implemented in different variations on 
different levels of sophistication in conjunction with various matrix-element generators or al- 
ready in full-fledged event generators, see e.g. [13,53-62]. Despite their differences they exhibit 
an assuring level of agreement [63]. 

In all these new approaches, parton emissions from matrix elements at a given perturbative 
order have to be balanced with corresponding emissions from a shower algorithm. Intuitively 
one may anticipate that dipole-like kinematics, leaving all particles of the splitting on their mass 
shells, may facilitate simpler procedures for this balancing. Furthermore, concerning matching, 
the CDM seems to be the more natural partner to the matrix-element part of calculations based 
on a subtraction method using antenna factorization [40]. 

In this publication, therefore, an extension of the "perturbative" dipole shower [16] as imple- 
mented in ARIADNE [8] to truly perturbative initial-state radiation is proposed, in contrast to 
the original ISR Lund CDM. Hence, the goal is to formulate the QCD evolution of a hard pro- 
cess initiated through a hadronic collision entirely perturbatively as a sequence of colour-dipole 
emissions. In particular, emissions associated to the initial state are treated as to directly emerge 
from colour dipoles spanned by the external parton lines. The beam remnants are kept com- 
pletely outside the perturbative evolution, their connection to the evolved cascade is left to the 
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hadronization to deal with. As a direct consequence, three types of dipoles and, hence, of asso- 
ciated radiation contribute to the full development of the final cascade, namely emissions from 
initial-initial (II), final-initial (FI), and final-final (FF) dipoles. Consequently, the emissions are 
denoted as initial-initial, final-initial- and final-state radiation (ISR, FISR, FSR), respectively. 
In order to model ISR and FISR in the fully perturbative version of the CDM proposed here, a 
backward evolution of the initial-state related radiation pattern of the shower is mandatory and 
automatically necessitates the inclusion of parton distribution functions (PDFs). 
The outline of the paper is as follows: in sec. [2] the basics of the dipole-shower model as im- 
plemented in ARIADNE will briefly be introduced. In addition, the treatment of ISR through 
final-state dipole splittings involving the beam remnants will be discussed. In the next section, 
sec. El the basic ideas of the newly proposed dipole shower are highlighted, especially its different 
ansatz for ISR simulation. This includes the generalization of the kinematical framework and of 
the evolution variables. In the following three sections, the kinematics and single-emission cross 
sections of all dipole splittings in various configurations of initial- and final-state partons will be 
detailed, recapitulating the case of dipoles consisting of final-state partons, the case implemented 
in ARIADNE, in sec. HI The kinematics and splitting functions of the new dipole types present 
in this model are discussed in sees. [5] and El The complete shower algorithm will be presented 
in sec. [7] and its performance for various physics processes will be highlighted in sec. [8l 



2 The Colour Dipole Model 

2.1 Physical background of the CDM 

The Lund colour dipole model (CDM) has strong connections to the semiclassical method of 
virtual quanta [64-66] , which equates the electromagnetic energy flux associated with the fields 
emitted by fast moving charges with an energy flux of equivalent photons. Owing to the large 
Lorentz boosts of the charged emitter, the corresponding electric and magnetic fields are or- 
thogonal to each other and they populate a plane orthogonal to the direction of motion of the 
emitter only. This amounts to a pulse of electromagnetic energy, given by 

dl(ui,b) ~ hadbJ \ — , (1) 

where b denotes the impact parameter, i.e. the distance w.r.t. the emitter; u is the frequency of 
the field component. It is bound from above by oj < p/mb, where p and m are the momentum 
and mass of the emitter, respectively. Equating this energy pulse I with a number of equivalent 
quanta n, 

dl = hujcLu, (2) 
and replacing the impact parameter with transverse momentum yields 
, a dk 2 , duj 

dn ~ — . (3) 

7T UJ 

A similar result emerges when considering bremsstrahlung off a charged particle, changing its 
otherwise straight direction of motion through a sudden "kick", or connected with the pair 
production of charged particles. Then, in the Breit-frame of the former process, or in the 
centre-of-mass frame of the latter, a rapidity y can be defined w.r.t. the axis of motion of the 



6 



2 THE COLOUR DIPOLE MODEL 



charged particle(s). A short calculation based on a full quantum-mechanical treatment shows 
that, neglecting spin effects, the number of bremsstrahlung-photons is well approximated by 

dn = — — dy , (4) 

VT U 

cf. [67]. Here, the rapidity must satisfy 

\y\ < \yo\, (5) 

and yo is the rapidity of the emitter (s). Rewriting energy through transverse momentum, 

k± = -V, (6) 
cosh y 

then leads to 

dn = - dy . (7) 

VT /Cj^ 

Because of its equivalence to eq. (j3J) , this equation exhibits the dominance of soft radiation in the 
semi-classical limit. In this context it is worth to note that the same limit is used in eikonal-type 
factorization of matrix elements employed, e.g. , in antenna subtraction methods [40,68] for the 
calculation of perturbative higher-order corrections to scattering cross sections in QCD. 
The simple formula for the semi-classical limit of photon radiation off a charged dipole, eq. (fT|), 
can be refined through a full quantum- mechanical treatment, including spin effects, see also later 
sections. However, the dominant features of the radiation pattern are already fixed by the simple 
formula, which in turn denotes the starting point for a shower simulation based on individual 
dipole emissions. The differential probability for such an emission to occur in an interval dp\ 
and dy is related to dV = da/dao given by 

dV^^^dy. (8) 
vr p z ± 

Here p± denotes a transverse momentum, which in the CDM is constructed out of Lorentz 
invariant quantities. Numbering the momenta of the particles after emission such that the 
newly emitted particle is labelled with "2", and, denoting the momenta before and after the 
emission with pi and pi, respectively, the emission can be symbolized as 

Pi + P3 = P1+P2+P3- (9) 
The squared invariant masses of sets of momenta are denoted as 

Sij... = (Pi+Pj + ■ ■ -) 2 and = (Pi + Pj + ■ ■ -) 2 ■ (10) 

A Lorentz invariant transverse momentum can be defined as 

p± = = — — , (11) 

S123 Sl3 

in agreement with [16,67] and the ARIADNE implementation. Moreover, a rapidity can then be 
computed through 
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7 



The Lorentz invariant choice guarantees a frame-independent description of the dipole splitting 
process. Using p± as the ordering parameter for subsequent emissions, a Sudakov form fac- 
tor encodes the non-emission probability between two scales p\ high and p\ low in analogy to 
conventional parton showers: 




In this form, the leading logarithms are resummed to all orders. 

The Sudakov form factor constitutes the basis of the simulation of parton emission also in the 
framework of the CDM. In contrast to ordinary parton showers, however, here the relevant 
objects are colour dipoles, which emerge naturally when considering the large Nq limit. In this 
limit, colour charges in the fundamental representation (quarks and antiquarks) have one colour 
partner, and colour charges in the adjoint representation (gluons) have two colour partners. The 
dipoles are built from pairs of such colour partners, and the emission of a gluon off a dipole 
effectively amounts to splitting the dipole into two. 

This self-similar process of dipole splitting, which is described in a probabilistic fashion, is 
easily encoded as a Markovian process in form of a computer program. Adding in the leading 
logarithmic behaviour and colour coherence as a dominant feature of QCD emissions results in 
a strict ordering of subsequent emissions such that the actual p± of a dipole splitting sets the 
maximal pj_ for the splittings of the two resulting dipoles. 

2.2 Initial-state radiation in the original CDM 

In ARIADNE, the only complete CDM implementation so far, initial-state radiation off incoming 
partons is not explicitly taken into account. Instead, ISR is redefined as FSR, where dipoles 
are spanned between potential final-state partons and the outgoing hadron remnants [24,69-71]. 
Considering DIS of leptons on hadrons, it can be argued that, as the hadron is in a bound 
state, all radiation originates from the colour dipole between the struck point-like quark and 
the hadron remnant - being an extended object composed of individual valence quarks and sea 
partons. Thus, an extended "antenna" is formed and from the electro-magnetic (semi-classical) 
analogy it follows that radiation of wavelengths smaller than the extension of the antenna is 
suppressed. Therefore, the original CDM was modified such that only a pj_-dependent fraction 
a(p±) of the remnant enters the splitting process [24]: 

«(^)=(£f, (14) 

where /i parametrizes the inverse size of the remnant and a refers to the dimensionality of the 
emitter, both being parameters to be tuned to data. In e + e~ annihilation the ("triangle") 
phase-space boundaries are approximated by \y\ < ln(M/p±), which now are supplemented by 
the extra condition y < ln(M/i/p^_). This obviously limits the range of accessible p± values in 
the splitting of the dipole of mass M. The strategy of sharing the recoil in such cases was inspired 
by the Lund string model, where an extra kink on the string (hadron) is interpreted as an extra 
gluon. This led to the introduction of recoil gluons to compensate for the recoil momentum 
associated with the part of the hadron, which participates in the emission. H Moreover, in cases 
where a sea quark is hit, the picture experiences further minor modifications. Taken together, 



2 For the ARIADNE rewrite in C++, the possibility of discarding recoil gluons completely is under consideration. 
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a good fraction of phenomenological, non-perturbative modelling enters the Lund CDM for ISR 
through all these assumptions. 

Next, consider Drell-Yan-like processes; there, a quark-antiquark pair annihilates to produce 
a lepton pair. In conventional parton showers, the two incoming quarks would emit secondary 
partons, typically simulated in a backward evolution algorithm [9, 12]. The recoil of these 
emissions is transferred to colour partners and the final-state leptons. In contrast, in ARIADNE 
the incoming quarks do not radiate but rather the two beam remnants, which are the only two 
coloured final-state objects before radiation (cf. left panel in fig. [T]). Then, the recoil of the first 
emission is compensated for by the final-state leptons [55,71], for all further dipole emissions, 
additional recoil gluons are added, if the emission occurred in phase space significantly away 
from the vector boson [71] . A further obvious refinement is the correction of the first emission 
to the corresponding matrix-element expression. The sharp phase-space cut-off is then replaced 
by a softer suppression function, in order to describe the high transverse-momentum spectrum 
of the vector boson. 

3 New approach to initial-state radiation using colour dipoles 

The principles underlying the proposal of this paper for the construction of a purely perturbative 
colour-dipole model are: 

• the maintenance of the probabilistic interpretation of emissions as encoded in the Sudakov 
form factor, which will be obtained from exponentiating single-emission differential cross 
sections; 

• the large Nq limit of the radiation pattern, and the restriction to account for the leading 
terms only, i.e. the leading dipoles, of this expansion; 

• the generalization of the kinematics and the evolution variables used in the original CDM 
and in ARIADNE to the case of ISR and FISR; 

• the factorization of the emission phase space and matrix elements around the soft limit 
(the radiation pattern has to be factorized in terms of 2 — > 3 splittings, to be derived for 
II and FI dipoles); 

• the utilization of crossing symmetries for the determination of dipole splitting functions; 

• the construction of on-shell kinematics for each splitting on an emission- by-emission basis, 
allowing stop and restart the cascading after any individual emission; 

• the backward evolution description of radiation related to incoming partons, and, conse- 
quently, the emergence of PDFs in the shower algorithm in a way similar to conventional 
parton showers. 

A number of issues are not at all covered here, which are, however, straightforward to include 
in some future work, namely 

• the comparison of different forms of splitting cross sections; 

• the approach's extension to account for finite, non-zero quark masses (here, all partons 
are treated as massless); 

• an extension to Supersymmetry; 



3.1 Kinematic framework 
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Figure 1: The Lund CDM approach to initial-state radiation in Drell-Yan processes (left 
panel) vs. the direct, perturbative approach as suggested by the new dipole-shower model 
(right panel). The treatment in modelling a first gluon emission is illustrated. 



• the QED radiation off the dipoles. 

To exemplify the impact of the principle outlined above, consider Drell-Yan processes; in contrast 
to the Lund approach, see sec. 12.21 i n the new dipole picture the primary dipole q~\q[ is directly 
formed by the two incoming quarks, and the emission will be calculated from the competition 
between gluon (see fig. [1]), quark and antiquark bremsstrahlung. The real-emission matrix- 
element information will directly enter, through the corresponding dipole splitting functions, 
as in the FF counterpart of emitting a first gluon off the qq' dipole. The boson's transverse 
momentum will be naturally generated, because the new initial-state momenta will be oriented 
along the beam direction. In case of an actual gluon emission, a system of two colour-connected 
successor FI dipoles emerges, namely a q\g^ and a gtq[ dipole. A further gluon radiated into the 
final state will then create a first FF dipole. If a quark is produced first, an FI dipole qtg\ and 
a successor II dipole g\q[ will be generated with again a dual role played by the gluon, here g\. 
Hence, the QCD evolution of the leading-order Drell-Yan pair production process eventually 
will involve all possible dipole types. 



3.1 Kinematic framework 

The occurrence of new dipoles and corresponding splittings is an immediate consequence of the 
suggested new CDM approach. A list summarizing the principal dipoles of QCD is shown in 
tab. [TJ Dipoles are labelled by k£, thus, in the most general way a splitting triggered by the 
emission of a (new) gluon g is expressed as 

kl -> kgl. (15) 

The notation is chosen such that the flavour and colour flow of all particles is outgoing. Three 
types of gluon emission emerge, related by crossing symmetry; any such splitting will leave the 
number of initial-state partons constant: 

gluon emission , 

quark emission, provided that k = q\, (16) 
antiquark emission, provided that 1 = q\. 

Here, the subscripts indicate whether the gluon emerges in the initial or final state. In the 
former case, this requires to replace the initial (anti)quark of the original dipole by the initial 
gluon and emit the corresponding antiquark (quark) in the final state. 



ki 




10 3 NEW APPROACH TO INITIAL-STATE RADIATION USING COLOUR DIPOLES 



QCDdipoles, kl 


II dipoles, i' i 


FI dipoles, / i 


FF dipoles, / /' 








mi 


<H9i 




9\9\ 


am 


gm 



Table 1: All dipole types appearing in QCD (the supplemental indices i or f label whether 
the parton is in the initial or final state, respectively; if clear from the context, the index f 
will be left out). The primes indicate that different quark flavours may constitute the dipole. 



Having clarified the notations used for the dipoles and their splittings, the kinematic objects will 
be introduced. First of all, the momenta are defined as incoming/outgoing if they are associated 
with the physical initial/final state. Those before and after the emission are denoted by p„ and 
p m , respectively, such that, expressed through the momenta alone the dipole splitting process 
can be written as 

P'k + StPi ► ft Pk + ft Pg + QPi- (17) 

Here and in the following the signature factors <fm = ±1 and q m = ±1 for partons in the final 
(+) and initial (— ) state. The before- and after-emission total momenta po and po then read 

-ftPo = ?- k Pk + ^Pl, (18) 

~% PO = ft Pk + ft Pg +S£Pl, (19) 

with the requirement that p$ = p\. Furthermore <fo = ft, an d the signature factor i.e. the 
association of the total momenta with the initial or final state is chosen such that the after- 
emission configuration refers to a production = —1, FF dipoles), scattering = —1, FI 
dipoles), or annihilation (^o = 1, II dipoles) process. Consequently, the four-vector po then 
corresponds to the four-momentum of the decaying parent dipole having mass \M\ such that 

Pi = M 2 = -Q 2 = p 2 , (20) 

with Q 2 arranged to be positive definite for FI dipole emissions. Accordingly, Lorentz-invariant 
energy fractions w.r.t. p$ are defined through!! 

x m = %^L. (21) 
Po 

The squared invariant masses of two- and three-parton systems are denoted by 

Smn — iSmPm ~\~ ^nPn) and Smnr — iSmPm ~\~ ^nPn ~\~ ft Pr) (22) 

where the inclusion of po and the expressions related to the momenta before the emission are 
understood. Concerning all gluon emissions considered here, the identity 

M 2 = s kg + s g£ + = s kg£ = -Q 2 , (23) 

holds true in general, since all partons are consistently treated as massless. 



3 The notion "energy fraction" is clear in the centre-of-mass frame of a parent FF dipole, where x m = 2E m /M. 
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3.2 Towards generalized evolution variables 

Next, the dipole evolution variables have to be generalized such that all emissions of all dipole 
types can be treated on equal footing and embedded in a consistent CDM-like evolution. The 
generalized variables should have the property of leaving the well-established FSR treatment 
unchanged and they should satisfy the constraint that all splitting cross sections, i.e. those 
involving initial-state partons as well, will follow the approximate form given in eq. ([8]). This 
would just manifest the universal features of QCD radiation in the soft limit, reproduced by 
eikonal distributions factorizing off the squared matrix elements in emissions off initial and final 
states alike, 

1 ( Pk Pi Y 2s kl ^ 



2 \PkPg PgPl) S kg Sgt 

Note that the right-hand side of this equation explicitly assumes massless partons. Following 
eq. (fTTD . the factor 2/p]_ b ecomes identical to the eikonal factor in the soft limit, and the collinear 
limits manifest themselves in the two-parton squared masses appearing in this p\ definition. 
Then, the generalized kinematic variables should exhibit the same singular behaviour in the 
soft /collinear limits for all dipole types, reflecting the crossing symmetry. Therefore, in this 
paper a generalized transverse momentum and rapidity are proposed in the form 



A 



$kg Sgl 



Skgi 



and 



y = - In 

y 2 



Skg 



(25) 



(26) 



Here, the invariant calculated including the signature factors <r m n r , i.e. through 

eqs. (|22|) . Clearly, for FF dipole cascading, all invariant masses are positive and hence the original 
CDM evolution variables of eqs. (jlip and (|12p are trivially recovered. For the other cases, the 
generalized form suggested here is a minimal Lorentz-invariant extension, guaranteeing a frame- 
independent evolution of the colour dipole. Moreover, these shower variables allow a global 
simultaneous ordering of all emissions. Given these generalized definitions, the identities 

\s k g\ = \M\p ± e- y and |vl = \M\p ± e +y (27) 

are found, indeed showing the similarity to the original Lund CDM. 

So, having introduced the kinematic framework, in the following sections, sees. HI [5] and El the 
derivation of the splitting kinematics for each dipole type (FF, II, and FI/IF) will be always 
pursued in four steps: 

• First, the evolution variables p± and y are identified. 

• Then, the limits of the emission phase space are deduced, which for the rapidity typically 
read y_ < y < y+. They guarantee that the evolution takes place within the physical 
region of phase space. These limitations are imposed through constraints on the evolution 
variables and, thus, determine the Sudakov form factor, see eq. (fTB"|) and sec. El 

• Together with the strict limits, approximate ones are stated, denoted e.g. for the rapidity 
by K_ < y < Y+, which allow an analytical evaluation of approximate Sudakov form 
factors. They are used in a Monte Carlo procedure, which finally corrects for the true 
form of the exact Sudakov form factors by means of a veto algorithm, see e.g. [3]. 
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• The on-shell three-parton kinematics of the splittings characterized by the central variables 
p± and y are constructed from the original two-parton configurations. Remaining degrees 
of freedom are fixed with a few additional assumptions, i.e. through four-momentum con- 
servation and a splitting-specific recoil strategy. The Lorentz invariant definition of the 
evolution variables guarantees the frame-independence of the actual construction. 



3.3 Dipole splitting cross sections and functions for QCD radiation 

In order to construct a parton shower as a Markovian process, the emission of any additional 
parton has to be factorized from the radiation of partons produced so far, such that the full 
radiation pattern can be built as a sequence of individual, mostly independent emissions. For 
the actual construction of a dipole shower, the individual parton emission should be modelled 
as being coherently shared between the two partons forming the dipole. The asymptotic form 
of these 2^3 dipole splitting cross sections has been presented in sec. El cf. eq. (JHJ), which 
constitutes the limiting case of soft emissions. In order to extrapolate to harder regions of 
emission phase space and to include spin effects, the differential cross section for a k£ dipole 
developing into a kgl colour-connected state is more conveniently written as 

dV M _ k l = = — D~ ki _> k t (p ± ,y) -fdy. (28) 

This defines the dipole splitting function D^i^^ in analogy to the splitting kernels employed 
in conventional parton showers. In both cases, the splitting kernels of dipole or parton showers 
incorporate refinements, which go beyond the corresponding eikonal or collinear approximation, 
respectively. Here, the D k i^ kg i may be deduced by analyzing differential cross sections for 
additional real emission of partons in comparison to the corresponding Born level processes. 
This then yields the single-dipole phase-space and matrix-element factorization, which has to 
work at least in the singular domains of the real-emission phase space. Accordingly, first-order 
real corrections are fully or partially encoded in the splitting kernels automatically. 
For the timelike case, the reasoning outlined above is realized by starting from the three-body 
decay rate of an object with mass M, 

dT 0^M> = \ M °^f9f'\ 2 d ^fgf'(P0; Pf,Pg,Pf>) , (29) 

where, in the massless limit, the dipole phase space and matrix element factorize according to 

dsfgds ii dip 

d^fgf (Po ; Pf, Pg, Pf ) = d%^ ff , (po =Po;Pf, Pf ) 16?r2 M 2 ( 3 °) 

and 

\ M 0^fgf'\ 2 - 8 ™sC D ff ,^ fgf , \M ^ ff , | 2 , (31) 

respectively. Here, taking Nq = 3 for the number of colours, the colour factor labelled C has 
been introduced explicitly. It typically takes one of the following values, C = Cf = = f f° r 

gluons emitted off quarks, C = Ca = Nq = 3 for gluons emitted off gluons, or C = Tr = | for 
gluon splittings into quarks. D denotes the dipole matrix element, which by definition correctly 
reproduces the singular terms of the parton emission process with potential differences in finite 
terms. Therefore, dT ^j- g j, can be expressed as 

Cot d(p 
dT o^f 9 f' - dr o^ff D ff'^fgf (P±. 2/> V>) d P 2 ± d y ^ (32) 



3.3 Dipole splitting cross sections and functions for QCD radiation 
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with px. and y taken from eqs. ([25]) and (J26|) . In most cases, the dependence of D on the 
azimuthal angle (p is neglected, thus, integrated out, such that the connection of the dipole 
splitting functions to the dipole matrix elements in the FF case reads 

D ff^fgf ip±,y) = iCp\b ff^fgf (p± , y) • (33) 

Note that here a gluon-sharing factor £ has been introduced because each gluon is contained by 
two dipoles. 

For the class of II dipoles, i.e. those consisting of colour-connected incoming partons, the ex- 
traction of dipole splitting cross sections has to be accomplished on the level of hadronic cross 
sections to correctly account for PDF effects and possible phase-space (suppression) factors. 
In this case, a 2 — > 2 scattering process rather than an 1 — ► 3 decay has to be considered. 
The differential cross section (using massless partons and having already integrated out the ip 
dependence) reads 

d<Tj/ i(g j)_, 0s (o g ) = h'( g ){x±,p F ) fi{x Tl fi F ) ^ \ M * <9*)-*09(0q)\ di dy cm , (34) 

where the usual 2^2 process Mandelstam variables, s = s^/j^i) and t = Sf g f qg .\ have been 
employed. Note that the quark emission case is signified by the parentheses. The fk(x±, p F ) are 
the PDFs. At leading order they can be interpreted as the probability of resolving a parton k 
inside the nucleon with light-cone momentum fraction x taken w.r.t. the nucleon's momentum; 
Pf names the factorization scale (defined in energy units), at which, pictorially speaking, the 
partonic substructure is probed. S and y cm denote the centre-of-mass energy and rapidity of the 
collider system, respectively. The 2 — ► 1 hadronic differential Born cross section for creating a 
particle of mass M through the matrix element Mi'iiqA^Q characterizes the t'i dipole's situation 
before the emission. Then, similarly to the FF case, the 2 — ► 2 matrix element squared can be 
cast into a factorized form, reading 

\-^t'i(gi)^0g(0q)\ 2 — 87rQ s C \-^t'i(qi)-*o\ 2 • (35) 

This hence allows to write the 2^2 differential scattering cross section in terms of the 2 — > 1 
Born term dov^^o where the tilde variables refer to the Born configuration: 



d&i'i(gi)^0g(0q) — ^fj/j^j)^ 



dy C m\ A'( s )(^±,Mf) fi{x T ,p F ) M 4 



dyemj fi'(q)(x±,PF) fi(x T ,flF) s 2 (p±,y) 

(36) 

Cct s * . . 2 

X t>T>i(q i i)^i>gi(qg i i){P±,y) dp ± dy . 

In contrast to the FF case, there is some additional freedom in arranging the actual recoils, since 
in principle the total energy of the splitting parton system will increase with each emission - 
additional momentum can be taken off the incoming nucleons. However, fixing the new centre-of- 
mass rapidity y cm removes this ambiguity. The choice in this paper is to ensure constant rapidity 
derivatives, thus, set up a recoil handling, which eventually shifts the original y cm through some 
function y that exclusively depends on the variables associated to the emission, 

Vera = y C m + y(M 2 ,s,i) . (37) 

Provided that these requirements can be satisfied, the II dipole splitting functions finally read 

n / \ fi'(g)(x±,VF)h(x T ,u.F) M 4 £Cpj ~ 

D mi i)^ g i( g9i i){P±,y) - /rw ( 5±|/ i F ) /i (fe F)AF ) p(p ±tV ) D m^'9i(mi)^y) (38) 
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and can be used to specify the associated differential splitting cross sections. In comparison 
to the gluon emission processes of FF dipoles, cf. eq. (j33[) . additional terms arise in each of 
the II dipole functions, namely a PDF weight, WpdFj which contains a ratio of PDFs taken 
at the respective momentum fractions and factorization scales before and after the emission, 
and a phase-space weight, Wpsp = M 4 /s 2 , which accounts for the altered incoming flux of the 
parton-level differential cross section. 

Similarly, the generic structure of final-initial dipole splitting cross sections dT > fi(fq i )-+fgi(fg i q) 
for gluon (antiquark) emissions into the final state is fixed through 

n , s fi(g){x±, fl F ) Q A jCp 2 L h 

Ji(q){x±,VF) [s(px,y) + Q 2 } 

where the Mandelstam variable s is related to the two-parton squared masses via s = Sf g (fq\- 
Note that the right-hand side of eq. (f3Uj) exhibits, as for emissions off II dipoles, the additional 
PDF and phase-space weights. The formula has been derived along the lines already employed 
for the II case, again in the limit of zero quark masses. As before, by comparing the differential 
hadronic cross sections before and after the emission, the emission part can be factored out when 
relying on dipole matrix-element factorization, 

IM)i(0gH/ s (/g)| 2 ^ 8na s C D fiifq . } ^ f gHfg.fi \Moi(Oq)^f\ 2 ■ (40) 



Like invariant squared amplitudes the dipole matrix elements D obey crossing symmetry. There- 
fore, as an alternative to the direct calculation of II/FI dipole matrix elements, cf. eqs. (|35h and 
(|40|) . these D expressions can be easily derived using crossing relations given that the FF dipole 
matrix elements have been worked out, see eq. (]3ip . In cases where there are different parti- 
cles in the final state, there is more than one possible crossing and, therefore, more than one 
corresponding dipole matrix element. The only remaining issues are the determination of the 
associated gluon-sharing and colour factors, £ and C, respectively. The latter are assigned ac- 
cording to the generic, large Nq, colour structure of the emission. 

Concerning dipole matrix-element factorization, there are two possible approaches to specify the 
dipole shower presented in this work. These approaches are: 

• Extract the D terms from the splitting cross sections employed in the Lund colour-dipole 
model [8, 15-17] for remnant-free dipole cascading. Results for II/FI dipoles are then 
derived from the corresponding FF dipole ones exploiting the crossing symmetry. In [16,72] 
the Lund differential cross sections have been shown to obey the correct QCD behaviour 
in the soft and/or collinear (Altarelli-Parisi) limit. El This reasoning therefore applies to 
the new cases as well. 

• Use the tree-level antenna functions presented in [41]. Their crossing symmetry has been 
exploited already while considering antenna subtraction with hadronic initial states, see 
[73]. Thus, the utilization of antenna functions, instead of the Lund kinematic functions, 
has the clear advantage of constructing a dipole shower out of subtraction terms that 
form the basis of the antenna subtraction method [40,41]; it therefore constitutes a very 
attractive alternative to the first approach. 

In the subsequent sections of this publication the first approach is being followed in order to 
allow for direct comparison with ARIADNE in the FF case. All relevant gluon emission types of 



4 Strictly speaking, this matching in the singular domains of QCD has been demonstrated omitting the influence 
of colour factors, i.e. it has been actually shown for the FF dipole matrix elements D. 
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2^3 dipole splitting functions will be listed and their parton-radiation characteristics in the 
various cases will be discussed. Finally, the performance of the full model is tested focussing on 
comparisons with experimental data. It is worth stressing, however, that the implementation of 
the second approach is straightforward and will be subject of a forthcoming study. 

4 Final-state colour dipoles 

In this section emissions emerging from FF dipoles are discussed. This is the traditional case 
already present within the original version of the CDM, implemented in ARIADNE. The dipole 
splitting process can be specified by 

f(k)f'(i) - f(k)gf'(i). (41) 

4.1 Final— final dipole single-emission phase space and kinematics 

Since the recoil of the emission will be completely shared between the three new partons, mo- 
mentum conservation, 

Po = Pf + Pj> = Pf+P g + Pf> = Po , (42) 

is realized between the momenta present before and after the emission. Note that apart from 
?o = —1, all other signature factors equal one. Neglecting parton masses, the relations 

< s mn = s 0r = M 2 {l-x r ) < M 2 , m ^n^re{f,gJ'}, (43) 

and the identity 

M 2 = Sf g + s g jt + s fj, , also expressed by 2 = Xf + x g + x j, (44) 

hold true. All energy fractions fall into the range < x r < 1, and, hence, the physics constraints 
imposed on the kinematic invariants s mn are satisfied. Following the steps outlined in sec. 13.21 
the {p\,y) phase-space parametrization can be characterized: 

• The two Lorentz invariant dipole evolution variables are 

P\ = Sj j^f - = M 2 (l-x f ,)(l-x f ), (45) 



cf. eqs. (fTTjl and (|25l) . and the associated rapidity y, 

. = ^n^ = ilni^, (46) 

2 Sfg 2 1 — X Jl 

cf. eqs. (fl~2j) and (f26|) . Therefore, the invariant masses can be re-expressed as, 

a gfl = M P± e +y , 
Sfg = Mp±e~ y , 

Sf p = M 2 — 2Mp± cosh y , (47) 



cf. eqs. (127]) . As expected, the dominant phase-space regions are characterized by p± — * 0, 
which points at p± , s utilization as the ordering variable. 
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The kinematic phase-space boundaries given through the relations in eqs. (|43p determine 
the (maximal) integration limits p\ high and y± stated in eq. (fT3|) . The determination of 
the precise boundaries is determined by the constraint 

Sfg + Sgj, = M 2 -s ff) < M 2 , (48) 

leading to the following symmetric rapidity limits 



M I M M 2 \ 

\y\ < arcosh = In \- \ — 5- — 1 • (49) 

11 " 2p ± \2 P± |4p2 j y J 

The largest possible value for p 2 ^ can also be read off these bounds, 
2 M 2 

P_L,max = —T- (50) 



Simple rapidity bounds overestimating the more exact interval are obtained, for example, 
s gf' 



from s f g , s 1, < M 2 ; this yields 



M M 

Y- = - In — < y < In — = Y+ , (51) 
P± P± 

which is nothing but the (y,z = ln^j) "triangle" commonly used to illustrate a dipole 
emission phase space. The effect of the sharper bounds now becomes apparent: they 
sizeably reduce the "triangle" area particularly in the central rapidity region, see fig. [2j 

Splitting kinematics: here the ideal frame to set up the new momenta is the centre- 
of-mass system of the parent FF dipole. Light-cone momenta!! w.r.t. the axis of the two 
original partons can conveniently be used. They yield 



p = (M, M, 0) - po = (M, M, 0) , 
p f = (M, 0, 0) Pf = (f±<*f,f±e-vt,f A 



0,M,0) - p f , = (/ie^/ieA/l) , 

Pg = Po-Pf- Pf ■ (52) 



~(') (') 

The f\_> and y f are specified by the particular recoil strategies that are used for the 
different types of FF dipoles. The choices taken here closely follow the approach presented 
within the Lund CDM, see e.g. [8]. Thus, for gluon emissions off qq dipoles, the Kleiss 
trick [74] has been implemented to treat the recoils: the (anti)quark will retain its direction 
after the emission with a probability x q(q)/( x q + x q)- For qg dipoles, the recoil of the 
emitted gluon will be compensated by the quark only. Specifying the kinematics of these 
cases (assuming, for example to preserve the direction of the /') leads to 

fx = ^^-sintf, f ± = (/ ± cos tp, fx sin ip) (53) 



5 In this work, light-cone momenta are defined as follows: q — (q+ ,q-,q±) where q± 
conditions can be intrinsically satisfied, if q — (mj_e z ,m±e~ z , q±) is chosen, using m\ = q 2 + q]_ and z = 
ln(q+/g_)/2 such that E q = m± coshz and q« — m± sinh^. 



4.2 Splitting functions for final-state QCD radiation 
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Figure 2: The phase space for gluon emission off FF dipoles; the dark-coloured region 
visualizes the available phase space. The bright colour is used to show the overestimation 
as given by the approximate limits w.r.t. the strict ones. Note that z — In ^j. 



and 

1 ,1 + costf / 0\ 

Vf = Q = 111 C0t 9 ' 54 

2 1 — cos ij \ 2 J 

where the polar angle i? is given through cos § = (x^ —x 2 f— x\ )/(2xfXp) and the azimuthal 
angle ip is taken to be uniformly distributed between and 2ir. Moreover, 

f± = Q, y'j = oo , preserving the product f'±e Vf = xpM . (55) 

For the distribution of recoils arising from gg dipole splittings, the simple specification will 
be corrected by rotating around the y axis in a way that the v\ of the parent gluons will 
be minimized, however, small perturbations introduced by an additional rotation around 
the x axis are allowed. 



4.2 Splitting functions for final-state QCD radiation 

In this section the refinements introduced by the Lund CDM [8, 15-17] to the simple eikonal 
splitting cross sections are briefly reviewed, in particular those for gluons arising from a qiqi 
dipole. Following the reasoning of sec. 13.31 the dipole splitting function for the 2 — ► 3 splitting 
qq' — > qgq' is worked out from the comparison of the real-emission process V — * qq'g to the Born 
contribution for the vector boson decay V — * qq'® For massless partons, the respective squared 
matrix elements averaged (summed) over colour and spin initial (final) states are 

"2~ ! 



\M r ^ g g\ 2 = 8 N c e 2 e\ p g Pg = 8 N c e 2 e\ 



N% — loo, x n "f" X n 



6 The inclusion of various correlations depends on exactly which processes are selected. For example, also the 
correlation of the leptons, producing the vector boson V, with the quarks could be accounted for by taking the 



processes ££' — > qq' (g) instead [74, 75] 
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The Lorentz invariant energy fractions of the emission are defined in eq. (|2ip . and, for simplicity, 
V = 7* has been chosen, see e.g. [76]. Therefore, 

dV^^ggg _ _ CfUs X 2 + Xg 

d Xq dx- q ~ 7 *^ 5 2vT (1-X g )(l-X f )' 1 j 

where = (JVq — l)/(2iVc) is the colour factor of this emission. Obviously, in this case the 
matrix-element factorization is exact allowing to read off the corresponding differential dipole 
splitting cross section according to eq. 



dV q ^ qgq = Cpo^ (I - ^ e+v) 2 + (I - ^ e^) 2 = Cpo, x\ (p± , y) + x\{p L , y) 
dp\dy 2n p 2 ± 2ir p 2 ^ 

The obvious overestimation of the exact result, 

J-papprox „ 

al qq^qgq _ <^FO<-s _f_ 

dp\dy 2tt p 2 ± ' 



(58) 



(59) 



in fact corresponds to a soft-gluon approximation neglecting quark spins, and allows for a direct 
implementation in a veto algorithm. Moreover, the dipole matrix element identified for the exact 
splitting can be cast into the following form: 

A _ x i "j~ x q _ 1 ( s i9 i s gq , ^ s qq s qgq \ 

U qq^qgq ~ 9 ~ _ I _ + + _ J ' V 0U J 

Vot b qg b gq J 

where the rightmost expression exactly reproduces the three-parton tree-level antenna function 
Ag(l g ,3 g ,2g) as stated in [41]. This case constitutes the easiest example for the compatibility 
of the two matrix-element factorization approaches. 

Similar reasoning can be applied to yield the splitting functions for gluon emission off quark- 
gluon and gluon-gluon FF dipoles [72]. It should be stressed, however, that in these cases 
the dipole matrix-element factorization is correctly achieved only in the singular limits of the 
emission. Taken together, the dipole splitting functions for gluon emission off final-state dipoles 
in the Lund CDM [15-17] read 



D ff<-+f 9 f> (p±,y) = f | Cr f | x n f f (p ± , y) + x™? {p± , y) 

- 2 k F A } c { F A } - D ?pT f9 p(p^y)- 



(61) 



They are all implemented in ARIADNE and will be used in the model presented here as well. The 
invariant energy fractions are given by 

x fJl = (62) 

Here and in the following, the parton-dependent exponents are defined as n qi9 = 2,3 and the 
curly-brackets notation is understood as | "' g°g e quark dipoles |. Note that £p = 1; in the splitting 

functions for dipoles consisting of at least one gluon the factor of £4 = \ enters, since gluons 
are shared among two dipoles. The D^ P p^ g p not only give upper bounds to the exact splitting 
functions, they also imply eikonal approximations to the splitting cross sections of eq. (|28p . 
A subtle issue in the formulation of a dipole shower is the assignment of colour factors. Obviously, 
for quark-quark and purely gluonic dipoles there are no problems, and, unambiguously, C = Cp 
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and C = Ca, respectively. For dipoles consisting of a(n) (anti)quark and a gluon, it is known that 
the colour factor cannot be pinpointed as straightforwardly as in the other cases, since e.g. for 
collinear radiation, the considered gluon emission can be traced back to either the (parent) quark 
or the (parent) gluon, such that in this limit the emission is therefore governed by Cp or Ca, 
respectively. Literally taken, these different colour-factor regimes have to be taken into account. 
This will lead to modifications of the corresponding dipole splitting functions and, possibly, to 
a decomposition (partitioning) of them into sub-contributions (subantennae) addressing these 
different regimes unambiguously. However, in the large Nq limit, underlying the construction of 
shower codes, this issue triggers sub leading effects only, since 2Cf, Ca — >• Nq keeping in mind 
that the gluon-sharing factor is £4 = |. 

5 Initial-state colour dipoles 

The first case, which goes beyond the original CDM, is radiation off an initial-state dipole ti of 
mass M. Two generic splittings based on gluon emission are available, namely 

i{k)i(l) -> l(k)gi{l) and %(k)i(i) -» q(k) gi i(£). (63) 
5.1 Single-emission kinematics 

Restating eqs. (fl~8l) and (fl~9l) for the II dipole scenario, i.e. setting % = q g i q = 1 and all other 
signature factors equal to —1, yields 

Pi'/qi+Pi = Po and Pr'/ gi +Pi = Po+Pg/ q , with pi = f Q = M 2 . (64) 

As already noted, the kinematics of the emission process here corresponds to that of a 2 — > 2 
scattering process rather than to that of a 1 — > 3 decay. The recoil of the emitted parton p g j q 
cannot be absorbed by p%'/ gi and Pi, since they are fixed to the beam axis. Thus, in contrast to 
the previously presented case, po 7^ po-lll For the scattering process, Mandelstam variables are 
defined as 

S = (PO+Pg/qf = {PV/g.+Pif = M 2 (l+X g/g ) > = S , 

t = (P0~p i ) 2 = (Pl>/ gi -Pg/q) 2 = M 2 (l-X t ) < 0, 

U = (P0-Pi'/ 9i ) 2 = (Pi-Pg/qf = M 2 {l-X V/g .) < 0, (65) 

where, again for massless partons, the bounds on s, t and u together with their parametrizations 
in terms of energy fractions are simple, and, furthermore, 

s + 1 + u = M 2 as well as / gi + X{ = 2 + x g / q . (66) 

As already indicated, the emission of a parton here requires an increase in s w.r.t. §q, related to 
an increase of the "Bj0rken-x". This is in contrast to the FF case, where the system's centre-of- 
mass energy remains constant. To deal with this issue, a generic parametrization is introduced, 
which relates the maximal partonic centre-of-mass squared energy to the squared mass of the 
parent dipole, 

Smax = aM 2 > s, such that 1 < s/M 2 < a < S/M 2 , (67) 
7 Unless the converse is enforced by a recoil handling in the 1'i dipole's rest frame. 



20 



5 INITIAL-STATE COLOUR DIPOLES 



where y/S is the centre-of-mass energy of the colliding hadrons. The limits on the invariants, 
detailed in eqs. (|65p and (|67p . clearly differ from the ones of the FF scenario, cf. sec. 14.11 This 
implies that the II dipole splittings arise in phase-space regions being distinct from the FF case 
and thus with a different kinematics. Consequently, the energy fractions populate new ranges 
compared to the FF splittings, viz. 

< x g/q < a - 1 , 

1 ^ x z'/gii x i — 1 "I" x g/q i 

2 < xy/ g . + Xi < 1 + a . (68) 

The phase-space parametrization is better worked out separately for both relevant II dipole 
splitting channels. 

5.1.1 Gluon emission phase space of initial initial dipoles 

First, the case of final-state gluon (g?) emission, i.e. t'i — > i'gi, is discussed: 

• The evolution variables are taken as suggested by eqs. (|25l) and (|26l) . The Lorentz invariant 
p\ thus reads 



A 



^ = M 2 (l-x,)(l-x ? ), (69) 



and the Lorentz invariant y is given by 
s 



y = — In 
y 2 



- In = - n , 70) 

2 t 2 1 - Xi ' K J 



such that the kinematic invariants can be re-written as 



s = 


SOg 


— 


= M + 2Mp_LCOshy 


i = 


SQi 


= Si'g 


= -Mp ± e~ v < 0, 


u = 




= S g i 


= -Mp ± e +y < 0. 



(71) 

The bounds on the Mandelstam variables - or equally well - on the invariant energy 
fractions translate, of course, into bounds on the evolution variables. As for FF dipoles 
emitting gluons, the more restrictive requirement is obtained from 

(a-l)M 2 = s max -M 2 > s-M 2 = -u-t = 2Mp ± coshy. (72) 

Hence, the allowed phase space, which is depicted in the left part of fig. [3j is described 
quantitatively through 

(a - 1)M . , 

\y\ < arcosh — , (73) 

2p± 



andd 

. (W-M 2 ) 2 (a-l) 2 M 2 



The maximal rapidity range is determined by the overall cut-off on p±: s max — M = 2Mpi jCUt cosh|j/ 



5.1 Single-emission kinematics 
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Figure 3: The accessible phase space for final-state (left panel) and initial-state (right 
panel) gluon emission off II dipoles. Bright colours indicate the phase-space fractions, which 
overestimate the respective allowed phase-space regions, which are shown in dark colours. 
The definitions of z are, z — In r^zjjM an< ^ z = m oM 9 { anc ^ 3' 1 emissions, respectively. 
Notice that the visualization of the g\ emission phase space is for a = 2. 



Weaker constraints are obtained from 



Smax-Af 2 > S-M 2 > -U,-t 

and, as in the FF case, they result in symmetric rapidity limits, 



\y\ < l n ' 



M 2 



In 



(a - 1)M 



Mp± pj_ 
These estimates again can be visualized by a "triangle" in the (y, z = In , p± 



(75) 
(76) 



plane. 



(o-l)Af - 

• The splitting kinematics will be detailed in sec. 15.1.31 together with that of the quark- 
emission process. 

Compared to the FF case, a new issue emerges: the maximal partonic centre-of-mass scale s max 
is not fixed and can be chosen. The actual choice then regulates the maximal size of the allowed 
emission phase space. This will be discussed together with the shower algorithm in sec. 

5.1.2 Quark emission phase space of initial initial dipoles 

Along the lines of the previous section, the phase-space parametrization and its consequences 
are now discussed for gluon emission into the initial state ((ft), i.e. (massless) quark emission into 
the final state: q\i — * qg\i. The details of the kinematics as outlined in sec. 13.21 are as follows. 

• The Lorentz invariant shower variables expressed through the Mandelstam variables read 



Pi 



V 



1 



In 



ts 

w 



- In-L 

2 -t 



M 2 ( Xi 



1 In 

2 x 



l+Xg 



1 



(77) 
(78) 



This allows to rewrite the kinematic invariants as 

+Mp ± e +y > M 2 , 
-Mp ±e - y < 0, 



s = 


SOq 


— s gi i — 


i = 


soi 


= S Q9i = 


u = 


SQ 9i 


— Sqi — 



M 



2Mp ± smhy < 0. 



(79) 
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implying, compared to the case of gf emission, a different shape of the valid (p±,y) phase 
space covered by this type of emission. 

The accessible rapidity range is (cf. right panel of fig. E|) 

, M [m 2 \ M aM 

In h a / + 1 = arsinh < y < In , (80) 

i o., ' 2p± p± 




where the left and right bounds result from u < and s < s max , cf. eq. (|67j) . respectively. 
This can be visualized in the (y, z = ln-^j^) plane as a "strip" that emerges in the point 
G/min, ^max = ~~ 2/min) and is confined between z = — y — In a and z = —y. The equations 
y m in = arsinh 5— ^ — and s max = Mp ±mSLX e y ™ in yield 

J_,max ' 

1 , 5 max 1 a , . 

ymm = 2 ln J—W = 2 ln ^I' (81) 

and 

Plmax = (W - M 2 ) ^ = o(a — 1) M 2 . (82) 

• The allowed phase-space region is safely covered by a "half-triangle" described through 
ymin < y < ~z. Accordingly, Ay = -z - y min = hi(p^ max /p5_)/2. 

• The splitting kinematics is presented in the next subsection. 

Finally, notice that, as for gi emissions, the single-emission phase-space maximally available is 
determined by the actual value given to s max . 

5.1.3 Construction of the splitting kinematics 

In the model proposed here the initial-initial dipole kinematics is directly constructed in the 
lab-frame. Particularly, to handle the recoils for the case of q[q\ dipoles, the strategy according 
to Kleiss [74, 75] has been implemented. 

Lab-frame kinematics: the fixed orientation of incoming partons implies that the emitted 
parton's recoil will directly be transferred to the entire final-state system, i.e. to all QCD and 
non-QCD final-state particles that are present before the emission takes place. As an example, 
consider the first emission in a Drell-Yan process, where the corresponding recoil is compensated 
for by the lepton pair. This recoil transfer results in po ^ po, and, therefore, a Lorentz trans- 
formation T defined through po = Tfio is necessary and will be applied on all particles (whose 
vectors are summed up in p~o). For the construction of the momenta, a light-cone decomposition 
w.r.t. the beam axis is well suited, such that, for massless partons, the situation before and after 
the emission is summarized as 



pi = yx + v S, 0, 0J — > pi = yx + v S, 0, 

Pi'/qi = (0, X-"JS, O) p V / g . = (0, X-VS, 0) , 



po = (M , M e~y° , 0) ->■ p = (m ± e yo , M ± e~ yo , -I 

P g/q = (t±e v <,£xe-*>,ty . (83) 



5.1 Single-emission kinematics 
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Furthermore, sq = x + X-S = M 2 and yo = y cm = ln(x + /x_)/2 with y cm denoting the centre-of- 
mass rapidity of the parton system. The x±, here functions of M, S and y cm , parametrize the 
momentum fractions of the partons w.r.t their respective hadron. Employing M"j_ = M 2 + i\, 
after the emission they read 

I e ±ye +M±e ±yo _ 
x± = -j= > x±. (84) 



Clearly, emissions leading to x± > 1 must be rejected. The vector £± = (£± cos ip, £± sin ip) and 
the quantity y e denote the transverse momentum and the rapidity of the emitted parton w.r.t. 
the beam axis, respectively. In terms of the Mandelstam variables, cf. eqs. (I65D . they are: 

?l = i and & = JLj-l-l-Z). (85) 

The azimuthal angle ip can in first approximation be assumed to be uniformly distributed, and 
s, t, u are determined by the evolution parameters p\ and y through eqs. ([71]) and eqs. ([79]) for 
gi and g\ emissions, respectively. The squared lab-frame transverse momenta are exemplified 
below as functions of pj_ and y. For gluon emission into the final state, 

f = M 2 pj = v\ 

± s 2 P± M- 1 coshy + 1 ' V ; 

whereas for quark emission, 

^ = M 2 p 2 L _ {M 2 + \t\)\t\ = {Me ~y ) 2 {2p±M -l sinhy _ 1) _ (87) 

When comparing both equations for the same ratio pj_/§, it becomes apparent that the emissions 

of quarks yield smaller lab-frame transverse momenta than those of gluons. 

To fix the last degree of freedom, an additional assumption is necessary, which is to preserve the 



rapidity of the system of outgoing particles, yo = Vo = Hem 
hand, s = x + X-S and y cm = ln(x + /x_)/2 = ln(n/t)/2 + y e 



^Having the complete emission at 
l®\ In more detail, 



1 , u I -t(s + u) \ 1 M 2 — t 
Vcm = Ucm + 7: m -x + In I - = y cm + - In — ; — j- , 

2 1 \^M 2 siu + {tu) 2 ) 1 s + t 

which exposes the impact of the yo = Vcm choice and shows that the system undergoes a rapidity 
shift during splitting. In addition, the new momentum fractions x± can be written down, 




(?) 

Finally the momenta, Pq , of all final-state particles, numbered by j, have to be transformed 
in order to account for the non-trivial change of po — » pq. Here, the Lorentz transformation 
T is specified as follows: the particles are boosted into the original dipole's centre-of-mass 
frame, afterwards the boost that forms po out of (M, ) is applied on them likewise. Altogether 
Pq = £>(— Po/pq)B(pq/pq) Pq = T p^ is computed. This finalizes the construction of the 
on-shell kinematics of an individual emission. 



If J/cm = ycm was naively exploited, the ratio of momentum fractions would remain constant, x+/x- — x+/x-, 
which constitutes a rather strange behaviour, since, for instance, very asymmetric starting configurations would 
persist to the end of the shower evolution. 

10 Particularly, for gluon emissions into the final state, y C m — ye = V- This simply expresses that rapidity 
differences are invariant under boosts along the beam axis. 
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5 INITIAL-STATE COLOUR DIPOLES 




Figure 4: Relevant Feynman diagrams contributing to vector boson production in associ- 
ation with a gluon (left panel, showing the t and u channel graphs) or with a quark (right 
panel, visualizing the i and s channel graphs). The modified colour flows due to the emission 
are illustrated by the thick lines. 

Improved description of lepton hadron correlations (Kleiss trick): when analyzing 
eqs. (|83|) again, it is noticed that, apart from the azimuthal angle 99, which eventually fixes 
the vector £±, all unknown variables are determined by Lorentz invariants plus the additional 
assumption yn = y C m-El]l n a first approximation, the choice is to uniformly distribute in azimuth 
w.r.t. the lab-frame, but more sophisticated schemes can be introduced correcting this simple 
ansatz. One such scheme can be derived from the work presented in [74] where it has been shown 
how to exactly factorize the first order tree-level corrections to the electroweak production of 
quarks. The corresponding Monte Carlo algorithm in fact is employed within the Lund CDM 
to arrange the splitting kinematics of q{q' { dipoles. In [75] this factorization was proven for 
scattering and annihilation processes involving initial states and corresponding algorithms were 
developed. Accordingly, for the q[q\ dipoles of this model, the suggestion of [75] has been 
employed to improve the splitting kinematics: the new momenta are constructed in the original 
dipole's rest frame in a distinct way, then they are transformed to the lab-frame such that the 
O-particle's rapidity is preserved. The essence is that the primitive ip choice is substituted by a 
prescription, which e.g. in Drell-Yan processes correctly accounts for correlations between the 
radiated parton and the leptons. As before, the particles associated to the parent II dipole have 
to be transformed, however they now undergo a more complicated series of transformations out 
of the (before-emission) lab-frame, i.e. more accurately 

# = *„ f A = P0 > + -/Z P0 - ) *W ( ^44 ) An f k) # , 

(90) 

where starting from the right, one applies to a momentum: the longitudinal boost into the 
dipole's rest frame, the alignment boost followed by the rotation that brings the newly incoming 
partons onto the light-cone axis maintaining the initial ± assignments, and the final longitudinal 
boost to satisfy that yo stays the same as it was before the emission, i.e. yo = yo- 

5.2 Initial— initial dipole splitting functions 

The first QCD-type emission in vector boson production (real-gluon bremsstrahlung or QCD 
Compton scattering) can be described as a coherent emission of a gluon or a(n) (anti)quark 
off the primary q^q[ dipole, cf. fig. [TJ For gluon bremsstrahlung qq' — ► Vg (<# emission) and 
QCD Compton scattering gq' —* Vq {g\ emission), the amplitudes can be worked out from the 

^Therefore, when neglecting the angle (p, it makes no difference whether the kinematics is arranged in the 
parent dipole's rest frame or in the lab-frame. 



5.2 Initial-initial dipole splitting functions 
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Feynman diagrams depicted in fig. 01 The partonic squared matrix elements, with colour and 
spin indices averaged (summed) over initial (final) states, can then be expressed in terms of the 
Born amplitude squared as 



and 



2iV c M 2 /2 fii 



LM^ p - I 4vra M 4 + u 2 — 2si 

respectively, whereas the Mandelstam variables have been defined in eqs. (|65|) . Again, in this 
particular case the factorization is exact and the equations above fix the dipole matrix elements 
Dq { q'.^q ig qi.(qg iq '.), which are in fact related to Dqqi_> qgq -t by crossing symmetry. The colour factor 
for the <7f emission is Cp. However, if the radiated gluon is assigned to the initial state, it 
actually is incoming and splits into a qq pair with one of the quarks entering the hard process. 
So, the colour averaging changes relative to the gi case by Nq/(Nq — 1); therefore, the colour 
factor amounts to Tr = |. 

For all initial-initial dipoles containing gluons, the dipole matrix elements may be obtained 
either directly in a similar way, or they are inferred from their final-final counterparts of the 
Lund CDM using crossing relations, D = cross D FF , cf. sees. 13.3 1 and 14.21 The recoil strategies 
presented in sec. 15.1.31 lead to trivial rapidity Jacobians, = 1. The initial-initial dipole 
splitting functions of eq. ([38]) are then fully specified: for gluons <?f emitted into the final state, 

n / \ h>{x±,VF) fi(x T ,u, F ) x¥'(p±,y) + x?(p±,y) 

L>T'i^T>gi{p±,y) = -pj-- — ~ \ f t~ — w F \ c l F \ 7 — ; s ', \ — 772 

f?{x±,p F )fi(x T ,u. F ) Ui UJ [xt( p± ,y) + x i (p ± ,y)-lY 

(93) 



where the energy fractions are given as 

xr,i(p±,y) = l + ^e^, (94) 
and, for gluons g\ radiated into the initial state, 



n , A f a (x±,y F )fi(x T ,fiF) „ x 2 q {p ± ,y) + x^{p ± ,y) 



fq(x±,M /<(5 t ,Af) [l + x q (p±,y)] 



2 



(95) 



where the energy fractions are characterized by 

x q ,i(p±,y) = Tl + j±e±y. (96) 

In both cases the overestimations _tPPP rox finally determine eikonal approximations to the im- 
proved splitting cross sections. The A/*pdf factors denote estimates for the respective upper 
bounds of the PDF ratios. The tilde variables refer to the before-emission state. 
All splitting functions discussed here are finite, i.e. the (soft and collinear) singularities of the 
various differential splitting cross sections defined through eq. (|28p are entirely contained in the 
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6 DIPOLES FROM FINAL-INITIAL COLOUR FLOWS 



l/p\ term of eq. (I28p . This nicely confirms that each eikonal cross section encodes the full sin- 
gularity structure of the exact result. For gf emissions, the invariant transverse momentum will 
tend to zero in either of the collinear limits that the gluon can have with the parent partons, i.e. 
the t or u variables turn independently to zero, or in the soft limit where x g — > and therefore 
t and u collectively approach the limit at zero. For g\ emissions, the divergence pattern is not 
as rich as for g$ emissions off II dipoles, since s is bounded to stay well above zero owing to the 
mass of the parent dipole. So, it only is critical if the emitted quark becomes soft or collinear 
with the incoming splitting gluon g[. No other radiating dipole contributes to this singularity, 
therefore £ = 1, consequently being omitted in the corresponding formulae above. 
The colour-factor assignment is unproblematic for quark dipoles q\q[ (see above) and also for 
gluon dipoles g\g\, where C = Ca and £, = £,a = 0.5. For II dipoles with a single gluon leg, the 
ambiguities beyond the large Nq limit appear in the same way as in the FF case. The following 
choices are currently made: final-state gluons are emitted adopting the Lund CDM choice of 
C = Ca (and £, = £,a = 0.5); for initial-state ones, C = Tr (and £ = 1) is selected adopting the 
result from the calculation for q\q[ dipoles. Since s > M 2 > 0, the selection C = T R at least 
ensures the correct behaviour in the singular limit t — > of <?i emissions. 



6 Dipoles from final— initial colour flows 

The branching of an FI dipole, fi, caused by a gluon may occur again in two ways by either 
radiating it to the final state, or to the initial state, releasing an antiquark instead: 

f(k) i{f) - f(k) g i(£) and f(k) <b(£) - f(k) 9i q(£) . (97) 



6.1 Single-emission kinematics 

Factorization implies that in deep inelastic scattering the evolution of the QCD particles proceeds 
completely independently of the evolution of the leptonic part. Therefore, not only the squared 
momentum transfer q 2 = —Q 2 from the lepton to the parton, probed by the scattering of the 
virtual photon, is a constant, but also q 11 = ke — k<t remains unaltered, when emitting QCD 
secondaries. This is used as the paradigm for the construction of the FI dipole kinematics in this 
model. Hence, similar to FF and in contrast to II dipoles, here the subsystem kinematically fully 
decouples from the rest of the cascade. Thus, in the FI case the partons directly participating 
in the splitting are affected only. Therefore, 

Po = Po (98) 

and 

Po+Pi/qi = Pf, Po+Pi/ 9i = Pf+Pg/q, with p 2 = M 2 = -Q 2 < 0, (99) 



such that Q may be interpreted as the "mass" of the parent dipole. Taking eqs. (I18p and (I19p . 
the signature factors are = = 1> all other ones equal —1. The underlying 2^2 

process implies to define the kinematic invariants for radiating FI dipoles as 

S = ijPO+Pi/g-f = (Pf+Pg/q) 2 = -Q 2 ( l + Xi/ 9i ) > = S , 
t = ijPO-Pff = {Pg/q-Pi/gf = -Q 2 (^-Xf) < 0, 

U = {pO-Pg/qf = (Pf-Pi/gf = -Q 2 (l-X g/ g) < 0, (100) 



6.1 Single-emission kinematics 
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where the identification of the energy fractions and the bounds are again given for massless 
partons. The Mandelstam variables then satisfy 

s + i + u + Q 2 = such that 2 + Xi/ g . = Xf + x g /q . (101) 

In analogy to the case of II dipoles, the maximal s is parametrized in terms of Q 2 as 

Smax = clQ implying that ^ s ^ s max < §. (102) 

Here, the quantity § = (po + P) 2 plays the role, which the squared collider energy S does for 
II dipoles, namely representing the maximal upper bound. The use of po = po and the rigorous 
definition of the Bj0rken-x variable, 

Q 2 



X B 



2p P 



where P labels the momentum of the incoming hadron, lead to 



Q [ 1 • 

xb 



(103) 



(104) 



This signifies that the Bj0rken-x determines the maximal range for the parameter a, namely 
< a < 1/xb — 1- Since parton masses are neglected, so = (Po + Pi/ qi ) 2 = and the Bj0rken-x is 
the momentum fraction x of the original incoming parton, Pi/ qi = x B P- Employing Pi/ gi = xP, 
it is found that xb < x = —Xi/g.XB < (a + 1) xb < 1 and the limits on x i / gi are clear: 



1-a < 



C i/9i 



1 + Xi/g. < Xf, X 



1 — a < Xf + x 



f t" Xg/q 



< -1 



< 1 



< 1. 



(105) 



The various other bounds then follow from eqs. (|10U|) and (|1U1|) . 



6.1.1 Gluon emission phase space of final initial dipoles 

First, FI dipole gluon emissions emerging into the final state, fi — ► fgi, are discussed according 
to the steps outlined in sec. | 



The evolution variables are identified as before by specifying eqs. (j25[) and (|26p for the 
case at hand. They read 



s fg s gi 



s fgi 



st 

^Q 2 



Q 2 (\ Xl \-l)(l-x f ) 



and 



2 



'9' 



s fg 



2 s 



2 Xi 



x f 



1 



(106) 



(107) 



whereas, using eqs. (|100|) . the Mandelstam variables can be rewritten as 



s 
i 
u 



soi = Sfg = +Qp±e y > 0, 
-Qp±e +y < 0, 



so/ 

so ff = Sfi = -Q 2 + 2Qp± sinh y < 0. 
Obviously, the rightmost relations for s and t are trivially fulfilled. 



(108) 
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1 


\ 




Az = ln(a+1) / 





Figure 5: Phase-space visualization of the final-state (left panel) and initial-state (right 
panel) gluon emissions off Fl dipoles. Again, dark colours indicate the available phase space 
for the emissions, whereas the fractions of phase space stemming from the overestimations 
are shown in bright colours. The definitions are, z = In ^ and z = In ^jphg for gt and gi 
emissions, respectively. Note that, for the visualization here, a — 1 has been assumed. 



• The largest phase space available is found from s < s max , cf. eq. (|102p . and u < 0, hence 

— In — — < y < arsinh — . (109) 

In the (y,z = ln^j) plane, see left part of fig. these bounds manifest themselves in a 
deformed "triangle", whose right side is curved to the inside diverging for y — > while 
approaching z = —y — In a for y — > do. The left side of the "triangle" is described by z < y 
and the intersection is at (y = z max , z m ax = In yl + l/a ), suggesting that 

Plmax = a(a + l)Q 2 = (s max + Q 2 ) 3 -^. (110) 

Similar to II dipole splittings, the maximum size of the emission phase space is dictated by 
the choice of s max , see eq. f|102j) . This can easily be understood, since the emission implies 
a new initial state with a larger momentum fraction taken off the corresponding hadron. 

• The exact rapidity interval is overestimated through the "triangle" bounds, which read 
z < y < -z + 2z max , resulting in Ay = ^{p± jiaax / P±) ■ 

• Again, the construction of the momenta is separately detailed, see sec. 16.1.31 



6.1.2 Antiquark emission phase space of final initial dipoles 



In this subsection, the phase-space parametrization for gluon emission into the initial state, i.e. 
antiquark emission into the final state, is discussed for massless (anti)quarks, fqi — > fg\q: 



According to eqs. (|25p . (I26|) and (jlOOH the evolution variables read 

S /3i S 9iQ 



and 



In 



S f9i 



ut 



ii4 

2 u 



Q\l-x q ){l-x } ) 



- In- 
2 1 



x f 



(111) 



(112) 
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The Mandelstam variables are then cast into the form 



s = 


s 0gi 


= s fg 


= 2Qp ± coshy - Q 2 > 0, 




i = 


so/ 


= S M 


= -Qp±e +y < 0, 




u = 


SQq 


= S /9i 


= -Qp±e' y < 0, 


(113) 



where the inequalities for t and u are satisfied by construction. 

The requirement < s < s max = aQ 2 in conjunction with eqs. (|113[) leads to 

arcosh^- < |,| < arcosh g - + Q2 = arcosh ( a + l )Q (114) 
2p ± ~ m ~ 2Qp ± 2p ± ' V ; 

where the inner and outer bounds follow from the lower and upper limits of the acces- 
sible s interval, respectively. Concerning the former the central rapidity region becomes 
unaccessible for emissions of p± < Q/2 (cf. right part of fig. [5]). In the latter case the 
similarities to gf emissions off II dipoles, cf. sec. 15.1-H are fairly obvious. If the available 
squared energy s max is fully used to generate the transverse momentum then coshy = 1 
and the maximal p± is achieved, given by 

P±,max = Jq 2 = j( a + l ) Q ■ ( 115 ) 

As in all cases involving initial-state partons, the adjustment of the size of the phase space 
is triggered by the choice of s max , which will be discussed in sec. [71 

Loose constraints stem from s max + Q 2 > —t, —ii and yield an increased phase space w.r.t. 
the precise one discussed above: 

\y\ < l n %±^! = l n ( Q + 1 )Q . (116) 

11 " Qp± p± v ; 

This invokes the usual "triangle" interpretation in the (y,z = In r^+h q ) plane. 
The splitting kinematics will be discussed in the next subsection. 



6.1.3 Construction of the emission momenta 

The basic construction principles mentioned throughout sees. [3] and 15.1.31 are, of course, taken 
over when explicitly establishing the FI splitting kinematics. The kinematic decoupling, encoded 
as po = po, alleviates the task, since Lorentz transformations will only affect the local splitting. 
Thus, including the fact that the squared dipole momentum p\ = (pf — Pi/ qi ) 2 = —Q 2 , the 
original dipole's Breit-frame constitutes a suitable frame to set up the three new four-momenta. 
Using light-cone notation again, in this frame the momenta read 

Po = (-Q, Q, o) 

p l/qi = (q, o, o) 

p f = (o, Q, o) 



Po 

Pi/9i 
Pf 

Pg/g 



-Q, Q, 0) , 

-x i/g . Q, 0, 



(l-x f )(-x i/g . - 1) x f 



\ x i/gi\ 
Po+Pi/ 9i ~Pf 



C i/9i 



\ x i/gi\ 



Q, b 



(117) 
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The Breit-frame transverse momentum is given through 6j_ = (b± cos if, b± smtp), where 
Vs i u f ~ ~ ~ Q 

b± = j^q2 = y( 1 - x f)(- x i/ 9i - 1 )( x f- x i/9 i - 1 ) ~\^~\- (118) 

Note that, for gt emissions, it becomes zero for u — > (in this limit the rapidity value associated 
to this emission coincides with the y+ bound, cf. sec. Ifi.l.ip . This just happens independently 
of the actual value for the evolution variable p±, therefore, ordering the emissions in p± does 
not impose any ordering in b±. Finally, the new Breit-frame momenta are transformed into the 
lab- frame. El 

Using eqs. (|117p the recoil strategy can directly be read off: before and after the splitting the 
initial-state parton is fixed to the + direction of the light-cone decomposition, therefore to the 
beam axis! 13 !, leaving the recoil to be completely compensated for by the final-state particle. 
Of course, more sophisticated recoil strategies following the ones of the Lund model and/or 
the Kleiss idea are possible, but not yet implemented. Especially the prescription for quark 
scattering processes given in [75] seems very attractive, since it includes the correlations between 
leptons and partons associated to the lowest-order DIS process and the first emission. 



6.2 Final-initial dipole splitting functions 

For qq[ dipoles emitting gluons, the two respective matrix-element factorizations of eq. (|40p 
can directly be specified utilizing the results^ for the two typical real-correction processes to 
leading order DIS, namely the QCD Compton and the boson-gluon fusion processes. These 
are compared to the sole scattering of a quark caused by a space-like vector boson; hence, the 
dipole's gluon emission will again be treated coherently and, moreover, exact factorization is 
achieved as in all other quark-quark dipole cases. This yields the corresponding dipole matrix 
elements, both of which in fact reflecting the crossing symmetry of the D qq /^ qgq / term: 

2 _L P - 9-r,n 2 x 2 , + X 2 n X 2 , + X 2 n 



_ 1 s* + t*-2uQ* _ x ^ x g _ x g[ + x g 

' Q 2 -st Q 2 {l + x q[ )(l-x q ) p\ - 1 :)) 



for <7f radiation off the qq[ dipole where C = Cf, and 



1 u2 + i2-2sQ 2 _ x^ + xl _ x 2 ,+x 2 q 

U ^m' Q2 Q2 {l _ x _ l){l _ Xq) p 2 . 



for (ft radiation off the qq[ dipole with C = Tr. For the definitions of the Mandelstam variables 
etc., see eqs. (|100p and the previous section. 

The dipole matrix elements of the FI dipoles containing gluon(s) are calculated either following 
the above procedure, or, alternatively, exploiting the crossing symmetry of the respective FF 
dipole matrix element taken from the Lund CDM. This completely determines the final-initial 
dipole splitting functions as introduced in eq. f|39|) : for emissions, 



n / x fi(-XiX B ,fj, F ) \x f (p±,y)\ n i +\xi{p ± ,y)\ n * 

fi{x B ,M t*/ x 2 {p ± ,y) 

<^ DF { { ; } C { ; } { 2(a 2 +1) } sD — M,V), 



(121) 



12 This is done by inverting the transformations that (1) align the lab-frame momenta pi/ q , and Pf and (2) 
rotate them afterwards onto the z axis. 

13 Note that the choice Pi/ 9i = —Xi/ gi pi/ qi = (—Xi/ gi XB\/~S, 0, 0) | lab f is Lorentz invariant. 

14 Colour and spin final-state summed plus initial-state averaged squared matrix elements are e.g. given in [76]. 
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with the energy fractions as functions of p± and y reading 

x u (p ± ,y) = ±l-?±e ± y, (122) 

and, for g\ emissions (for antiquarks emitted into the final state), 

n / x f gi (-x gi x B ,HF) „ \xf{p±,y)\ nf +x 2 g( P± ,y) 

B/.-Wpx.») = /(xg|W ^ 

(123) 



where the energy fractions are then given in terms of the evolution variables by 

X f,q 



(P±,y) = l ~^ e±V i such that x gi (p±,y) = -^^coshy. (124) 

The modulus ensures that the terms in the rightmost numerator of the exact splittings are pos- 
itive definite. Additionally, the eikonal approximations are displayed, which again overestimate 
the true form of the splitting functions. 

As in the previous cases, the dipole splitting functions are finite, such that the divergences are 
fully encapsulated in the l/p\ term. For gluons emitted into the final state, collinear/soft limits 
(t — ► or/and s — > 0) appear as before, where the (collinear) singularities for gluons are again 
only fully accounted through the inclusion of the contributions of the neighbouring dipoles (cf. 
the choice of £, £ = £a = 0.5). 

If the gluon is radiated into the initial state, the incoming gluon may split collinearly and, there- 
fore, in singular domains w.r.t. both the emitted antiquark (t — > 0) and the "other" final-state 
parton associated to the emission (u — ► 0), cf. eqs. (|100p . This is in contrast to the situation 
of II dipoles where a collinear divergence cannot emerge between the incoming gluon and the 
"other" parton, since in this case it belongs to the initial state. Turning to the discussion of the 
soft infrared limit, the gluon g\ itself cannot become soft, since it is coupled to the initial state. 
Therefore, t and u cannot vanish at the same time, i.e. the soft limit is kinematically shielded, 
which is also clear from eqs. (]105p . In case a soft antiquark is being emitted, a singular effect 
only occurs once it is also collinear with the splitting gluon such that t — * (the associated 
disappearance of s is non-singular). 

The colour factors are chosen similarly to the previous cases, with the same reservations concern- 
ing the collinear limits. For FI dipoles, not only a final-state gluon emerging from a quark-gluon 
dipole gives rise to the ambiguities, in this case, also the antiquark emission into the final state 
induces them on the same level. This is related to the fact that this splitting, as already men- 
tioned, is singular when either the emitted antiquark (t — > 0) or the "other" parton (u — ► 0) 
in the final state become collinear with the initial-state gluon, cf. eqs. (jlOOp . The ambiguity 
here occurs when this "other" final-state parton is a gluon, apparently resulting in a collinear 
splitting governed by Ca rather than Tr. 



7 The complete shower algorithm 

In this section the dipole-shower algorithm is presented, which models the full QCD radiation 
picture in terms of initial-state, final-initial and final-state colour-dipoles on purely perturbative 
grounds. This formulation of the shower aims at resumming effects at leading logarithmic 
accuracy while producing exclusive final states of partons. These are generated in a Markovian 
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process, iterating individual emissions. In analogy to conventional parton showers, a Sudakov 
form factor constitutes the central probabilistic quantity that determines the full development 
of the cascade. This will be discussed first, before the procedure of evaluating the evolution 
variables that characterize a single emission briefly will be explained. Finally, the showering 
algorithm will be fixed by specifying its relevant parameters and scale choices. 



7.1 The Sudakov form factor 

The evolution variables are given by the invariant transverse momentum p± and the invari- 
ant rapidity y, defined in eqs. ([25]) and (|26|) . respectively. Since the dipole splitting functions 
D(p±,y), cf. eq. (|28p. are finite throughout, the entire singular structure of each emission cross 
section in each case is incorporated as a term \jp\- The concept of "time" therefore is realized 
through p±_ , which thus operates as the (leading) variable ordering the emissions within the cas- 
cade. Consequently, y is considered as the associated variable. Within the generic framework, 
cf. sec. EJ all emissions are treated on equal footing, resulting in a competition between different 
available channels at each evolution step. The Sudakov form factor is obtained from integrating 
the corresponding differential single-emission cross sections d ^T dy (which are positive definite) 

in suitable boundaries of p± and y. Summing over all allowed splitting channels {kl — ► kg£} 
and exponentiating the negative result finally yields the Sudakov form factor: 



A (pi,st tJ pi) = exp< 



. PI 
p± 



where 



y+ (p±,a) 

164) = £ / iyn-^M.v)- P*> 

{M^kgl} y _(p ±) a) 

In this form the Sudakov form factor resums the leading logarithms as encoded in the dipole 
splittings to all orders, and, hence, can be interpreted as a no-branching probability. Accordingly, 
the two infrared divergent contributions of virtual and unresolvable real emission cancel each 
other below the infrared cut-off leaving an overall finite result. Thus, A(p^_ stt ,p±) quantifies how 
likely a state consisting of a number of dipoles will not emit any further resolvable parton between 
the start scale p\ stt and a lower (cut-off) scale p\. The quantity /xr denotes the renormalization 
scale in energy units for the evaluation of the (running) strong coupling. Typically /jr is given 
as a simple function of the evolution variables to include some higher-order virtual contributions 
beyond the leading-logarithmic approximation [77,78]. If reduced to the case of FF dipole 
evolution only, the expression for the Sudakov form factor of course becomes equivalent to 
eq. ()13p of the Lund CDM. Note that in the more general case, the rapidity limits y± also depend 
on the scaling quantity a, cf. eqs. (|67p and (|102p . i.e. on the choice of the maximal available phase 
space. Additionally, their actual functional form depends on the particular emission channel. 
For notational brevity, this has been omitted, but is clear in view of eq. (|49p . eqs. (|73p . (|80p 
and (I109p . (|114p . The presence of the PDF ratios in the I/FI splitting kernels naturally yields 
a Sudakov form factor including these ratios. This resembles the typical backward evolution 
treatment, where the ratio of parton densities ensures that the parton composition of the hadron 
is properly reflected in each evolution step [9] . 



7.2 Generation of the emission's Sudakov variables 
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Finally, the actual differential probability (the probability density) for some branching to occur 
at p\ then reads 

dP _ dA( P l stvP l) _ J(p 2 ) 

dpi ~ dpi ~ pi A ^tt,Pj- (127) 

Subsequent emissions are ordered in p±, i.e. their start scale pi stt is identical to the pi of the 
last parton radiation. This generates the Markov chain. Note that this still leaves the initial 
starting scale - dubbed initializing scale - for the very first emission, pi ini , to be selected. The 
choices made here are detailed in sec. 17.31 



7.2 Generation of the emission's Sudakov variables 

In the model a valid pair of evolution variables is generated by exploiting the strict p± ordering, 
which enables to treat any dipole and each of its emission channels separately. Therefore, for 
each single channel, a trial (pl,y) pair is generated according to its probability density 



dpi 2-ir pi 

y- (p± , a) 



J dyD~ kUkgi (p ± ,y) 



x exp < 



p i,stt y + (p±,a) 
dp L a s [n R (p±)] 



pi 2ir 

p 2 , ~~ y- (p± , a) 



(128) 

dyD~ kl ^ kg£ {p ± ,y) 



A valid {pl,y) pair generated according to the distribution eq. (|127|) is finally obtained by 
iterating over all channels picking the one of largest pi from the ensemble of all trial pi's. 
The procedure of selecting such a trial (pi , y) pair for a single dipole emission channel follows 
the standard Monte Carlo technique (hit-or-miss method) of the veto algorithm [3] exploiting 
that, for any given pair, the eikonal approximations gathered throughout sees. 14.21 15.21 and 16.21 
overshoot the respective dipole splitting cross sections. Each of which yields a fully integrable 
and invertible probability density, for which the (j>l,y) selection can be solved analytically 
using two random numbers. The respective simpler density can be easily read off eq. (|128j) 
when replacing a s [fi^], y± and D k ^ kg£ by a sufficiently larger a™ ax , loose rapidity bounds 
Y± (overestimating the actual rapidity interval) and approximate splitting functions D~ pT ° x , 
respectively. The correction to the true form of the single-channel density is then achieved by 
accepting the trial pair with a probability equal to the ratio (correction weight) reading 

w _ »s(wO D U-.kgi Ay(p ± ,a) 

afax D ^AY( P± ,a)- [LM) 

However, there are additional kinematical constraints, such as the demand for valid momentum 
fractions x± < 1 in II dipole evolution, see sec. 15.1.31 once violated, they translate into rejection 
of the trial emission, implying the generation of a new trial emission for the considered channel, 
starting over from the rejected p± value. 



15 Both forms are positive definite and describe differential cross sections. Therefore, employing them as kernels 
in Sudakov exponentials will always yield Sudakov form factors smaller than one, such that these form factors 
can be interpreted as all-orders expressions in leading logarithmic accuracy for emitting no parton between two 
evolution scales. 
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Note that the PDF-ratio overestimations, A/pDF; present in the approximate I/FI splitting func- 
tions are taken from a dynamically self-adapting table in order to improve the generation effi- 
ciency. The last term of above equation exhibits the correction for the exact rapidity interval, 
where Ay = y + — y_ and Al" = Y + — Y- . 

7.3 Scale choices, starting conditions and iteration principles 

Finally, the remaining free scale choices are fixed. This completely defines the (default) cascade- 
generating algorithm of this shower model. 



Renormalization scales: 



Mr 



,2 



FF = T ^ 



„2 J2 



= 2z(l-z)p 2 ± = P \ , r < ^ (130) 

I/FI v ' y± l + cosh(2y) ~ 2 K ' 



are used for the argument of the running strong coupling, where z is defined as a fraction of 
squared two-particle masses of the partons partaking in the emission: 

Kl _ 1 (131) 



\s k9 \ + \s g e\ l + e 2 y' 
An offset of 0(1 GeV) ensures the evolution to proceed well above the Landau pole Aqcd- 

Factorization scales: for a new (trial) emission ftp is calculated according to 

ft = (4fci) d / 2 /x F 2 - d , (132) 

where d = 1 and d = 2 are employed for II and FI dipoles, respectively. The modified transverse 
momentum squared, k"j_, is computed from the emission's p\ and y and the mass of the parent 
dipole. Based on eq. (|86p Ar^ is frame-independently defined as 

k 2 = \skgSgt\ = \M\p 2 ± 

1 \Skgt\ + \Skg\ + 1 8 gt \ \M\ + 2 P± COsh y ' 

and intended to function as a more natural scale for the argument of the parton densities, 
since it better compares to the lab-frame squared transverse momentum. The respective old 
factorization scales associated with the state before the emission are encoded in the p-p values, 
whereas, for the very initial case, /iF,ini is adopted from the hard process. 

Initializing scales: the subsequent cascading off the core process starts at the hardest scale, 
p\ = p^inj, which can not be set independently of the underlying process. Generally it should 
guarantee that the shower strictly evolves in the soft and collinear phase-space regions only. 
Here, the following choices are made for three different scenarios of hard 2 — ► 2 processes: 

• Showering off a single qq dipole as in e + e~ — > qq processes: the start scale is set by the 



squared mass of the parent dipole, p ± 



qq-prod. ^0 



M 



2 



Showering off a single q\q[ dipole as in Drell-Yan processes: here, a px,ini estimate is gained 
by inverting eq. (|133|) for y = 0, which yields 



P±,h 



DY 




(134) 
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and gives p_i_ ; im | DY = (1 + V%)M, provided that k\ max = §o = M 2 . For the same reason 
as in the previous item, this estimate is found by restricting the true transverse momentum 
of final-state gluons radiated off q\q[ dipoles, cf. eq. ([86]) . For vector boson production, 
in this model, the first emission is matrix-element corrected per construction. Therefore, 
the restricted scale may be discarded, the shower instead evolve freely with the initializing 
scale set as largely as kinematically allowed. 

Showering off a multi-dipole state as in pure QCD jet production: recalling the definition 
P± = \skg s g e\/\skgi\, all possible combinations for this fraction can be calculated using 
the strong particles provided by the hard process. The combination yielding the lowest 
p\ should represent a sufficient estimate for the initializing scale. Applied to QCD jet 
production, the minimal numerator is given by mm{ut, st, sii} employing the Mandelstam 
variables of the 2^2 QCD core process. For the denominators, a mean squared mass 
determined by |(s + \t\ + \u\) is used replacing the vanishing \sk g e\- Taken together, 



P_L,ini 



9 min{'ui, st, su} 

= ; — • (135) 

Jets s + It + it 



The assignment of large Nq colour flows is straightforward and unique for the first two examples. 
In the latter case, the most likely flow among the possible ones for a given set of partons will be 
determined by the actual kinematical configuration of the QCD 2^2 core process. 



Maximal phase space: for II and FI dipole regimes, the limits on the evolution variables 
vary with the choices of the s max parameters restricting the phase space of a single emission for 
decreasing s max values, cf. eqs. (|67|) and (|102H . The default settings impose no extra restrictions, 
they hence are s max I n = S (see sec. 15 . 1 j) and s max | FI = § (see sec. 16. ip and allow full access 
to the centre-of-mass energy as given by the collider. 



Cascading: each chain (colour-singlet), once appeared, is independently evolved, with the only 
potential exception of recoil transfer from an II splitting. This does not spoil the further evolution 
of the corrected chain owing to the Lorentz invariance of the shower formulation. Starting off 
P± ini consecutive emissions are decreasingly ordered in p 2 ^ within a chain. The dipole splitting 
that generated the largest p\ in a certain evolution step is finalized in its kinematics, its p\ is 
used as the new start scale for the trial emissions of the next round. The procedure continues 
until the infrared cut-off has been reached. 



Cut-off and hadronization aspects: the cut-off is always taken on p\, hence denoted by 
Pxcuf With a prescription avoiding the Landau pole in a s , it can, in principle, be chosen 
arbitrarily small, since then the setting of the renormalization scale is safe and the Sudakov 
suppression quenches the appearing soft and collinear divergences. 

After the cascading is finished, the interface to the hadronization - currently described through 
phenomenological models only - does not require any special treatment inside the shower. The 
conversion of the shower partons into primary hadrons proceeds similarly to the case of conven- 
tional parton cascades. 



8 First results 

In this section, the newly developed dipole shower is validated by comparing its predictions of 
QCD dynamics to data and other Monte Carlo calculations. To this end, the following physics 
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processes are studied: 

• the production of vector bosons and their subsequent hadronic decays in e + e~ collisions 
at LEP1 energies, 

• the inclusive production of Drell-Yan e + e~ pairs at Tevatron and future LHC energies, 
and, 

• the inclusive QCD production of jets at Tevatron energies. 

The shower model presented here has been implemented into the event generator SHERPA, and 
supplemented by an interface to the Lund string fragmentation routines of PYTHIA 6.2 [3], which 
are provided by the SHERPA framework. Cascading starts off the corresponding hard 2 — > 2 
processes, which are generated inside SHERPA utilizing its facilities of evaluating matrix elements. 
Only light-quark flavours, i.e. massless quarks, are considered. For the simulation of hadronic 
collisions, all predictions have been obtained from the CTEQ6L set of PDFs [79]. In accordance 
with the choice in the PDF, the strong coupling constant has been fixed by a s (Mz) = 0.118 
and its running is taken at the two-loop level. The dipole-shower cut-offs related to final-initial 
and initial-initial dipole evolution are both set to 1 GeV, i.e. p\ cut FI = p\ cut n = 1 GeV 2 . 
In contrast, p\ cut | FF is tuned by hand together with the Lund string model parameters. 
The lower panel in each of the plots presented below visualizes the (MC— reference)/reference 
ratio, where the "reference" (ref) is given by the data as long as they are available. The bright 
band always illustrates the uncertainty of the respective measurement. 



8.1 Hadron production in electron— positron collisions 

The testbed to exclusively validate the performance of the sole final-state piece of the dipole- 
shower modelEfl is the process e + e~ — > Z°/j* — > hadrons, where the qq pair produced in the hard 
process will initiate the cascade. The QCD Monte Carlo predictions can be compared with large 
sets of data, which, for example, are available from the LEP1 measurements. The data precisely 
test the QCD dynamics of hadronic final states produced at the Z° pole. The parameters of 
the shower and the hadronization model were tuned by hand, i.e. the value taken for the strong 
coupling at Mz was specified, the FF cut-off p\ cut | FF of the dipole shower was adjusted, 
and, suitable values for the Lund string model parameters a (PARP(41)), b (PARP(42)) and 
a q (PARP(21)) were found. The method employed for that is sufficient to yield first significant 
results. However, it can not be compared to the effort of delicate Monte Carlo tuning procedures 
as presented in [81] and foreseen in [82] in order to automatize the procedure. The "naively" 
tuned parameters read: 

a s (M z ) = 0.1254, 



P_L,cut 



FF 
a 



0.54 GeV 2 , 



0.29, ( 136 ) 
b = 0.76 GeV" 2 , 
a q = 0.36 GeV . 

Since massive quarks are not handled yet, the dipole shower always started off massless qq 
pairs. At S = M| a mean parton multiplicity of (A/" pa rton) = 9.24 and a mean charged-particle 



16 Note that, for pure final-state cascading, the splitting of gluons into quark-antiquark pairs already has been 
included, and the implementation is nearly identical to the treatment proposed in the Lund CDM. 
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Figure 6: Durham differential jet rates as a function of the jet-resolution parameter y cut ; 
dipole-shower prediction vs. DELPHI data taken from [80]. Here, the light-coloured band 
represents the sum of the statistical and systematic errors. 

multiplicity of (M c h) = 20.47 are found, where the latter is somewhat below the experimentally 
detected value of (A/" c h) = 20.92 ± 0.24 [22]. Figs. [6] and [7] show a selection of distributions 
obtained with the dipole shower and compared to DELPHI data taken at = 91.2 GeV during 
the LEP1 run. In fig. El Durham differential jet rates, Y n _, n+ i, are presented up to Y^q. They 
disentangle at which y cut = 2m.in{Ef , Ej}(l — cos 9ij)/S values an n + 1 jet event is merged 
into an n jet event according to the Durham jet clustering scheme [83]. The agreement with 
the data taken by the DELPHI experiment [80] is very good, in particular the description around 
the peak positions. All predictions tend to be somewhat below the bin means for low and high 
values of the jet-resolution parameter y cut . 

Event shape variables probe the pattern of QCD radiation for both soft and hard emissions 
arising from the primary qq dipole. Therefore, in fig. the shape distributions of 1— thrust, 
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Figure 7: The dipole-shower predictions for the event shapes 1— thrust, sphericity, thrust 
major and thrust minor vs. DELPHI data [22]. Again, the light-coloured band represents the 
sum of the statistical and systematic errors. 



1 — T, thrust major, T ma j or , and thrust niiuor, T m [ nor are displayed together with the sphericity, 
S, placed in the top right corner of the figure. The former are all obtained from a linear 
momentum tensor, whereas the latter stems from a quadratic one, thus, puts more emphasis 
on high momenta. All dipole-shower results are compared, once again, to DELPHI data [22]. 
The low-value parts, which are sensitive to soft emissions, are all quite well described, except 
for larger deviations in thrust major and minor. Differences in these observables also appear, 
even somewhat larger, for instance for the new shower of HERWIG++ [27] and the new shower 
presented in [36] based on Catani-Seymour dipole factorization. For a very recent comparison, 
please cf . [84] , where the value of the strong coupling constant has been determined at the Z° 
pole using the results from a first NNLO calculation for e + e~ 3 jets [85]. Although the 
soft parts of these distributions are all affected by hadronization corrections and their careful 



8.2 Inclusive production of Drell-Yan lepton pairs at hadron colliders 
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Figure 8: Boson transverse-momentum distribution in e + e _ + X as predicted by the 
dipole shower for two different choices concerning the initializing scale. The Monte Carlo 
calculations are compared with CDF data [86] taken during Run I at the Fermilab Tevatron. 
The right panel depicts the very soft region of the distribution only. 



modelling and tuning, the good behaviour of the dipole shower in describing soft emissions can 
be seen as a consequence of exponentiating the eikonal rather than the collinear limit of QCD 
radiation. The predictions for hard emissions agree somewhat worse with the data. The last 
two bins of the 1— thrust distribution are overestimated signalling a slight excess of spherical 
events, whereas thrust minor is underestimated for high values. 

Taken together, the agreement with data is satisfactory. This allows to conclude that the final- 
state piece of the dipole shower sufficiently is under good control. 

8.2 Inclusive production of Drell— Yan lepton pairs at hadron colliders 

In the scope of hadronic collisions, the processes pp(pp) — ► Z /^* — ► e + e~ constitute the simplest 
and cleanest testbed for the further validation of the dipole shower as they form the initial-initial 
dipole counterpart of the qq timelike evolution. 

Tevatron Run I predictions 

The transverse-momentum distribution of the lepton pair is heavily influenced by additional 
QCD radiation arising in both soft and hard phase-space domains. Owing to its clear signal, 
this spectrum has been measured with high precision by the Tevatron experiments. It is shown in 
fig.[8]for lepton-pair invariant masses in the range 66 GeV < M ee < 116 GeV. Two hadron-level 
predictions produced by the dipole shower are confronted with data from a CDF measurement 
at y/S = 1.8 TeV [86] and normalized to the experimental inclusive cross section. They differ in 
their choice of the initializing scale, using, first, p_i_ 5 i n i = {l+V%)M ee and, second, pj_,i n i = P±,max 
(cf. eqs. (JTIJ) and ([82]) ). In the latter case the shower evolves totally unconstrained, exploiting 
the fact that the first emission is corrected for the true matrix element by construction and may 
hence appear at a scale exceeding M e 2 c . This in turn sets the highest scale for all consecutive 
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DY production eV+X : boson rapidity 




Figure 9: Rapidity spectrum of the vector bosons in inclusive Z /^* production at Teva- 
tron Run II energies of = 1-96 TeV as predicted by SHERPA and the dipole shower in 
comparison to recent D0 data [87]. The grey solid, the black dashed, the green dashed and 
the red solid lines give the SHERPA CKKW ume = 1, CKKW n M E = 2, the default and the 
unconstrained dipole-shower predictions, respectively. 



emissions. The whole treatment eventually leads to a good agreement with the data for large 
pr- In contrast, the dipole shower with restricted initializing scale gradually undershoots data 
above px = 60 GeV before it dies off rapidly around pt = 80 GeV. 

The figure's right part contains a close-up of the peak region on a linear scale, almost identically 
predicted by both dipole-shower variants. The turn-on of the distribution is well described. 
Around the peak, narrower described by the data, a slight excess is found, followed by an 
underestimation of the data for the region above 12 GeV. The predictions include an additional 
improvement for very low pr, namely an intrinsic transverse-momentum smearing, which has 
been tuned by hand to these low pt data according to a Gaussian to have a mean (width) of 
0.3(0.4) GeV.|£0 Without this correction the shower p? spectra would slightly shift to the left. 



The assignment of an intrinsic transverse momentum to the hard process is a non-perturbative correction 
applied after the shower phase, and, therefore, irrelevant for shower kinematics. 



8.2 Inclusive production of Drell-Yan lepton pairs at hadron colliders 
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Figure 10: Left panel: impact of standard scale variations on the shower evolution ex- 
emplified by means of the pT,e distribution. Here, the blue lines reflect the uncertainty of 
the prediction from the unconstrained shower, in both panels given by the red solid line. 
Right panel: impact of the Kleiss correction on the transverse-momentum distribution of 
the e~. Here, the blue dashed and the grey dotted lines denote the outcome of the uncon- 
strained dipole shower neglecting Kleiss corrections and of the SttERPA CKKW time = 1 
merging, respectively. Both panels depict shower-level (SL) results (lacking hadronization 
corrections) . 



Tevatron Run II predictions 

One more validation against data is presented by considering the rapidity distribution of the 
decaying vector boson, where the Drell-Yan lepton-pair mass has been restricted to the interval 
71 GeV < M ee < 111 GeV. The QCD NNLO theoretical prediction for this inclusive observable 
has been calculated in [88] and very good agreement with data from a recent D0 measurement 
[87] has been observed over the full rapidity range. 

Here, hadron-level predictions are presented that have been obtained from the unconstrained 
dipole shower - denoted by "Dipole shower, max" in the plots from now - as well as from the 
P± ini restricted dipole shower, which is taken as the default, since the matrix-element correction 
of the first emissions does not apply beyond Drell-Yan processes. The comparison also shows 
SHERPA outcomes resulting from the CKKW merging of parton showers and tree-level matrix 
elements up to time extra partons. This merging method has been validated in many other 
comparative studies [53,63,89,90] or even to data [91]. Here, two such inclusive samples, for 
n ME = 1 and nME = 2, were generated with SHERPA using its version 1.0.10. All results are 
displayed in fig. [9] and confronted with the D0 data [87]. There hardly are any shape differences 
visible between the various Monte Carlo predictions. This nicely confirms that the II dipole 
kinematics is eventually well fixed by preserving the rapidity of the final-state particles, cf. 
sec. 15.1.31 However, compared to data, all Monte Carlo shapes are somewhat wider showing an 
excess of up to 20% for large rapidities. 

A rough estimate for the uncertainty of the shower predictions can be gained from varying the 
values taken for the and /iR scales within the shower algorithm. To this end, their defaults 
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were multiplied/divided by 2. The /xp scale enters through the PDF weight, and /Ur as the 
scale of the running strong coupling in the single-emission probabilities, cf. sec. 17.21 This is 
exemplified in the left plot of fig. [10] for the px e distribution, where the uncertainty band for the 
unconstrained (intrinsically first-order matrix-element corrected) dipole shower is shown to cover 
the same-order-of-accuracy prediction stemming from CKKW Z° + 1 jet merging. Therefore, 
both descriptions are in good agreement. The right part of fig. [TO] is the verification for the 
importance of the Kleiss corrections for emissions off q[qi dipoles. Their application yields a 
hardening and, therefore, an improvement of the single-lepton px, e spectrum of about 20%. 
SHERPA CKKW time = 1 again serves as a good reference, since it accounts for the full first- 
order lepton-hadron correlations. 

LHC predictions 

The correct energy extrapolation of the dipole shower is now verified by comparing various ap- 
proaches at LHC centre-of-mass energies. Therefore, the unconstrained dipole shower is studied 
w.r.t. SHERPA's CKKW merging for iime = 1,2,3 and SHERPA's pure showering realized by 
APACIC++ [14,92], which is a virtuality-ordered parton shower in the traditional sense resum- 
ming large logarithms in the collinear rather than the soft limit of QCD radiation. 
Most of the observables presented here require the exclusive definition of jets, which has been 
attained according to the Run II kx algorithm [93, 94] using the parameter D = 1 and an un- 
constrained T) range in order to include forward-jet effects. The jet px threshold is given by 
Prjet > 20 GeV. All distributions are simulated at the hadron level and normalized to unit area, 
which allows for direct shape comparisons. Many plots show CKKW predictions for hme > 1, 
which helps estimate the impact of describing the next-to-first extra parton emission by matrix 
elements as well. 

The pseudo-rapidity and px distributions of the e + e~ pairs are shown in fig.llll They are largely 
determined by the pattern of QCD emissions, in particular by the hardest one, and measure the 
recoil of the lepton pair against all other final-state particles. Hence, these inclusive observables 
are first defined beyond the leading order Drell-Yan pair production process. For the r\ ee spec- 
trum shown in the left panel, the maxima and the central rapidity region are somewhat more 
and less pronounced by the dipole shower than by the CKKW predictions, respectively. As it 
can be seen, the improved description of second-order emissions results in a further enhancement 
of vector bosons that are central in r\ space. In contrast, the APACIC++ prediction features a 
considerably larger dip in the central pseudo-rapidity region as a consequence of lacking suffi- 
ciently hard emissions, since this shower's start scale is constrained by the mass of the lepton 
pair to ensure evolving in the collinear and soft phase-space regions only. The right part of fig.fTTI 
contains the PT.ee distribution on a double-logarithmic scale to provide good insight to both soft 
and hard px domains. In the hard tail the dipole-shower result is 30% below the CKKW refer- 
ence; the difference in the low px part amounts up to 40%, whereas the dipole shower clearly 
puts emphasis on the soft region and predicts a slightly lower peak position. The agreement 
is still satisfactory and the deviations can be traced back, for the very soft part, to different 
parameter settings for the fragmentation of the partons (including intrinsic kx smearing), for 
the range 1 GeV < px < 20 GeV, to different radiation patterns generated by the dipole shower 
and the vetoed APACIC++ shower used for the CKKW merging, and, for the high px tail, to 
differences in choosing and processing the scales in both approaches. Typically, the inclusion of 
next-to-one extra parton emissions at the matrix-element level leads to a further increase of the 
high px tail. This is found in fig. [TT1 where the effects can become as large as 40%. 
Fig. [12] presents the jet pseudo-rapidity, r/j, and transverse-momentum, pr,ii distributions of the 
first three jets. These observables directly probe the jet structure of the events. For rji, the 
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Figure 11: Pseudo-rapidity spectrum (left) and px distribution of the lepton pair in inclu- 
sive Z°/7*+jets production at the LHC. The comparison is at the hadron level between the 
unconstrained dipole shower (green solid lines) and various SHERPA results, namely CKKW 
time = 1 (grey solid lines), taken as the reference curve, CKKW time = 2 (black dashed 
lines), CKKW time = 3 (red dotted lines) and APACIC++ (blue dot-dashed lines). 



dipole-shower predictions are quite similar and in all cases narrower w.r.t. the SHERPA CKKW 
predictions. For Pt%, the predictions of the unconstrained dipole evolution agree quite well with 
the respective ones of SHERPA CKKW time = 1> again on a 20%-40% level, confirming that the 
scale setting by the first (the unconstrained) emission reasonably constrains the subsequent one. 
The pt hardness of the jets predicted by the inclusive two- and three-jet merging is of course out 
of reach for the dipole shower. Such higher-order corrections can only be included by matrix- 
element parton-shower merging techniques or a matching with respective NLO calculations. 
The top left plot in fig. [13] depicts the vector boson rapidity spectrum obtained under the addi- 
tional requirement that the first and the second jet appear well separated in rapidity according 
to yi2/2 < — 2. Except for APACIC++ predicting a strong tendency of the boson to accompany 
one of the jets, all codes give flat spectra for central rapidities, and, remarkably, the dipole- 
shower result agrees well with that of SHERPA CKKW ume = 2. A similar pattern is found in 
the |Ary eej i| = \r/ ee — r)i\ distribution shown in the top right plot of fig. [T31 The dipole-shower 
curve hardly deviates from the CKKW curves, which reliably describe this observable owing to 
their higher-order contributions. This tellingly highlights the effects of the improvements of the 
dipole splitting functions beyond the eikonal approximation. In contrast, APACIC++ describes a 
suppression for low |A?7 eej i| values. 

In contrast to the shape comparisons above, tab. [2] provides some insight concerning inclusive 
and exclusive jet rates normalized to the total inclusive cross section. The results given by the 
unconstrained dipole shower are close to those of SHERPA CKKW jt-me = 1- Moreover, it is also 
found that the default dipole shower predicts much more soft jets compared to APACIC++. 
The large phase space available for additional QCD radiation at the LHC will lead to the copious 
production of jets. Here, this testbed provides an excellent means to study the effects of the 
s mM reduction (see sec. I7.3|) . In tab. [2] and the bottom row of fig. Q3] various predictions are 
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Figure 12: Pseudo-rapidity (left column) and transverse- momentum (right column) distri- 
butions of the first three jets in inclusive lepton-pair production simulated for LHC energies. 
Dipole-shower results are shown in comparison to those obtained by the CKKW method of 
SHERPA. Labelling is as in fig. QT] 
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Figure 13: The rapidity of the lepton pair for the first two jets satisfying j/ij/2 < —2 (top 
left) and the modulus of the pseudo-rapidity difference between the vector boson and the 
leading jet (top right). Both of which is simulated for inclusive Drell-Yan e + e~ production 
at the LHC. The dipole-shower outcome is compared to those received from various SttERPA 
runs. Labelling is as introduced in fig. 111! 

The bottom plots exemplify the impact of the choice of s max on the y ee for yiy2 < —2 
and cos #13 distributions as predicted by the unconstrained (solid green curves) , the default 
(dashed blue curves) and the dipole shower where s max is set dynamically (dotted red curves). 
The cosine of the angle between the first and the third jet, cos #13, is determined in the rest 
frame of the first-plus-second jet system. 
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0.009 
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Table 2: Cross section ratios as obtained from the various Monte Carlo approaches for 
inclusive and exclusive (hadron- and shower- level "SL") jet rates at LHC energies. Jets are 
defined according to the Run II kr algorithm [93,94] and required to have prjet > 20 GeV. 

confronted with each other, namely those of the unconstrained dipole shower, the default shower 
where p\ ini = (1 + V2) 2 M 2 e and the dipole shower denoted by "dyn" where additionally s max 
is set dynamically according to p 2 ^ max = p 2 ^ ini = (1 + \/2) 2 M 2 e (cf. e.g. eqs. flH]) and ([52]) in 
sees. l5.l7Tl and l5,1.2l respectively). The default shower loses hard emissions, such that the ratios 
of tab. [5] are somewhat smaller w.r.t. those of the unconstrained shower. For the "dyn" variant, 
the normalized cross sections decrease further, however, dramatically fall if a rapidity separation 
for the first two jets is required and imposed by yiy2 < —2. Furthermore, while the spectra 
presented at the bottom of fig. [13] show only mild differences between the unconstrained and the 
default shower, the predictions of the "dyn" shower deviate considerably: the y ee spectrum for 
2/12/2 < — 2 is far too peaked in the central y region, which also contradicts the performance of 
the CKKW references shown in the top left corner of the same figure. And, in the cos #13 distri- 
bution of fig. [13] the third jet is significantly less collinear w.r.t. the first and second jet, where 
the angle #13 between the first and third jet is taken in the rest frame of the combined first and 
second jet. Hence, in all examples, the s max reduction manifests in a suppression of forward and 
larger separated emissions (jets), which can be understood, since, firstly, the p\ max act as kine- 
matic upper bounds to all start scales p\ stt , in particular initializing scales p\ ini . Secondly, for 
reduced s max , the generation of large \y\ is suppressed for a single emission, eventually avoiding 
the production of sets of two-parton invariant masses, \s^ g \ = \M\p±e~ y and \s g e\ = \M\p±e +y , 
where one of them is very small and the other very large. Taken these findings together - also 
recalling the good agreement with the CKKW results - it is evident that using s max n = S and 
Smax I fi = § is a good choice. 

8.3 Inclusive jet production at hadron colliders 

The copious QCD production of jets is a typical and large phenomenon at hadron colliders, how- 
ever, from a theoretical point of view, the task of calculating and/or simulating these processes 
at higher orders in the strong coupling is more complicated and rather involved. Clearly, QCD 
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jet production severely tests the entire shower algorithm and goes beyond the tasks handled by 
the dipole shower so far. There are several reasons for this: the primary state is now given as 
a multi-dipole configuration formed by the 2^2 hard QCD processes according to their (large 
Nq) colour connections, including those that link initial- and final-state partons. Possibly, all 
dipoles form only one colour singlet or even a "gluonic ring" . Matrix-element corrections for the 
first extra emission in jet production are absent in the dipole splitting functions; this in turn 
requires to carefully constrain the initializing scale, such that the shower evolution is guaranteed 
to proceed in the soft and collinear limits of QCD emission only. 

To validate the dipole shower the observables listed below have been considered in more detail. 



Dijet azimuthal decorrelations at Tevatron Run II energies 

The dijet-decorrelation observable measured in the transverse plane between the two hardest 
jets, A0dij et = |0i — 02 1 j provides good insight to the occurrence of additional soft and hard 
radiation. There is no necessity to reconstruct further jets. The clear full-correlation signature 
given by A0dij et = it washes out in the presence of extra emissions. Since, the strength of the 
decorrelation rises in dependence on their hardness, the dijet decorrelation can be used to verify 
any candidate choice of setting the initializing scale. El 

The observable was subject of a recent measurement by D0 at Tevatron Run II with the data 
taken in different px.i = PT,m&x windows of the leading jet [95]. The details of the analysis are: 

• Reconstruct cone jets for R = 0.7, 

• require pr,2 > 40 GeV, and, 

• require central jet rapidities, 1 2/1,2 1 < 0.5. 

Fig. [Tl] shows the data overlaid with predictions for various choices of the initializing scale: be- 
sides the default given in eq. (|135p . two alternatives have been implemented, namely a geometric 
mean reading 



n 2 



2 . 2 S t U 



3 ^QCD = 3 , 2 , 1 2 , „, 2 » ( 137 ) 



Jets ~ ^ U " s 2 + t 2 + « 2 ' 

and a more enhanced scale defined as 



P_L,ini 



= (1 + V2) 2 ^ , (138) 
Jets S 



using the Mandelstam variables of the core process. The scale denoted by "II sc." is taken 
according to the latter equation, eq. (|138|) : this one denoted by "QCD sc." corresponds to 
eq. (|137l) . and that marked as "low default sc." follows from dividing eq. (|135l) by a factor of 
3. Obviously, the dipole shower initiated through the low default scale does not account for 
enough hard emissions and overshoots for soft ones. The other predictions are quite similar, 
with the "II sc." and "QCD sc." variants giving slightly harder and softer results w.r.t. the 
default case, respectively. The default performs best and will therefore be employed in all what 
follows. Its predictions still tend to undershoot the data around A^dijet = 2.8 in all PT,max 
windows of the leading jet, however, keeping in mind that some gluon splitting processes have 
not been fully taken into account yet, the agreement is satisfactory giving evidence that also 
other model-intrinsic scales, such as jlp, and /iR, have been chosen reasonably. 



18 Showers preferably should predict the distribution for small decreases of A^dijot = T, the tail may be corrected 
by the inclusion of matrix elements beyond order o? s . 
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Figure 14: The dijet azimuthal decorrelations in different pr.max ranges. Dipole-shower 
results for different choices of the initializing scale are overlaid by data taken by D0 during 
Tevatron Run II [95]. 



Dijet mass spectrum at Tevatron Run II energies 

With the p± t mi finding in hand, the dipole-shower prediction for the dijet mass spectrum is 
confronted with data measured during Run I by the D0 collaboration [96]. The analysis requires: 

• the reconstruction of jets using a cone algorithm with R = 0.7, 



• jet transverse energies above 30 GeV, and, 

• the dijet candidates to satisfy |t7x,2 I < 1-0- 
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Figure 15: The dijet mass spectrum as measured by D0 during Tevatron Run I [96] 
compared with the prediction of the dipole shower initialized at scales according to eq. (|135|) 
and for no phase-space restriction, i.e. using the default s max settings. 



As it can be read off fig. [T5| the comparison versus data with the dipole-shower results being 
normalized to the cross section observed in the experiment shows encouraging agreement. 

Test of colour coherence at Tevatron Run I energies 

An interesting measurement and analysis was carried out by the CDF collaboration during 
Tevatron Run I, searching for evidence for colour coherence mpp collisions at yS = 1.8 TeV [97]. 
Discriminatory observables were found for three-jet events featuring a hard leading jet and a 
rather soft third jet. They were shown to be sensitive to the correct treatment of QCD colour 
coherence in parton shower simulations. Here, similarly to the treatment in [36], this CDF 
analysis is used to test whether the proposed dipole shower is capable of describing the colour- 
coherence effects seen in the data.Lj The requirements of the CDF study read: 

19 Evolving in terms of colour dipoles is said to automatically account for soft colour coherence owing to the 
eikonal structure of the dipole splitting cross sections, however, the colour factor ambiguities for quark-gluon 
dipoles (discussed in sec. 14.21 ff ) require more serious investigation in this direction. The comparison with the 
CDF data is just a first step. 
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Figure 16: Colour-coherence tests in inclusive three-jet production at Tevatron Run I 
energies according to a CDF study presented in [97] : (left panel) pseudo-rapidity distribution 
of the third jet and (right panel) the angle a (defined in the text). Experimental errors are 
statistical only and the histograms are normalized to their respective binwidth. For the 
latter three observables, dipole-shower (shower-level) predictions under full (blue solid lines) 
and restricted (black dashed lines) emission phase space are shown in comparison with the 
(detector- level) data of the CDF measurement [97] . 



• Jets are defined through a cone algorithm, using R = 0.7, 

• the two leading jets are constrained to \r)\^\ < 0.7, 

• they have to be oriented back-to-back within 20 degrees, i.e. \<p\ — <j>%\ > 2.79, 

• jet Et thresholds have to be respected for the first jet and all next-to-first jets of 110 GeV 
and 10 GeV, respectively, and, 

• for the a angle only, a cut on AR23 = 1/ (772 — %) 2 + {<j>2 — fa) 2 nas to be imposed, namely 
1.1 < Ai? 23 < vr. 

• The angle a is defined through 

tana = ~ 3?) , (139 ) 

\fa ~ fa\ 

In fig. [16] the comparison between detector-level data and dipole-shower predictions obtained at 
the shower level is shown for the 773 and angle a distributions. As pointed out in [97], these two 
observables receive small detector corrections only, which is not the case for the AR23 separation 
of the second and third hardest jet in (77, eft) space. The latter is known to be strongly affected 
by detector effects, therefore, not considered here. 1^1 If colour-coherence effects are modelled 
correctly, 773 should arise broader and feature a significant dip for central values. The a spectrum 



For the same reason, in the study of [36] Ai?23 has not been taken into account either. 
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Figure 17: Pseudo-rapidity distribution of the third jet and spatial separation between 
the second and third jet in exclusive three-jet final states simulated for Tevatron Run II 
energies. The plots are obtained requiring a larger jet pr threshold for the hardest jet to 
pass. The SHERPA shower prediction generated by APACIC++ in CKKW scale-scheme mode 
(black dashed curves) is compared to the dipole-shower predictions, namely for default (blue 
solid curves) and lowered initializing scales (red dashed curves). 



should be minimal for small \a\ followed by a clear rise towards larger positive angles. As can be 
seen in fig. [161 the dipole shower predicts these characteristics, providing fairly good evidence 
that colour-coherence effects are reasonably modelled. The agreement with data deteriorates 
once the prediction is taken from a dipole shower where the s max setting has been (considerably) 
reduced. This again emphasizes that the natural choice is to assign the full phase space to single 
emissions by using the default s max settings. 



Exclusive three-jet final-state challenge 

Recent CDF measurements have found an excess in data of exclusive three-jet events with 
small Ai?23, which is not described by available tools, such as PYTHIA (Tune A) [98,99]. In a 
first qualitative study the potential of the new dipole shower to predict AR23 differently w.r.t. 
traditional leading-log showers is estimated. Therefore, the following analysis has been applied: 



• Require jet reconstruction according to the cone jet algorithm, use R = 0.4, 

• use general cuts on jets of pT,i > 20 GeV and \rji\ < 2.5, 

• additionally, use \r]i\ < 1.0 for the hardest jet, and, 

• consider the trigger-jet effect, i.e. demand px,i > 40 GeV. 

The results are presented in fig. [T7] and show the pseudo-rapidity spectrum of the third jet 
together with the Ai?23 distribution mentioned before. Concerning the former, the dipole-shower 
variants are found to generate steeper spectra, while inspecting the latter, the shower starting 
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from the lower initializing scale prefers populating the region of small jet separations the most. 
The APACIC++ prediction rises towards smaller separations as well but stays below the dipolc 
shower curves and features a broader tail. Additionally, by relaxing the pt,i > 40 requirement, 
treating all jets likewise, the higher px t i threshold was identified as a major source in projecting 
out the peak for low AR23. The different behaviour of the dipole and parton shower eventually 
can be seen as a consequence of generating the full radiation pattern differently. For example, 
recalling that the p\ definition includes a product of two-parton squared masses S12S23, a small 
s\2 with a potentially small angle between parton 1 and 2 can still be compensated by a large 
S23 giving the same p\ . In a 1 — > 2 splitting usually there is no such freedom of compensating a 
small S12, it might be rather cut away by the parton-shower cut-off. In conclusion, the tendency 
of the dipole shower to enhance the production of spatially less separated jets should be studied 
in more detail and more realistically including underlying event simulation etc. and, possibly, a 
direct comparison to data. 



9 Conclusions 



In this publication, the colour-dipole shower approach based on the Lund Colour Dipole Model 
[8, 16] has been extended. In this model, the description of sequences of QCD emissions is for- 
mulated in an expansion around their soft limit. In the context of hadronic collisions, a novel, 
perturbative description of initial-state showering based on the emission properties of colour 
dipoles has been developed, which is in clear contrast to the corresponding Lund ansatz. In 
summary, initial-state radiation is treated directly and not redefined by final-state radiation 
arising from final-state dipoles that contain extended colour sources, which, therefore, are sub- 
ject of a semi-classically motivated suppression of high p± emissions. In contrast, in the new 
model the hadron remnants are completely kept outside the evolution. The fully perturbative 
treatment led to the introduction of new dipole types, which contain incoming partons. 
The description of gluon emissions off colour dipoles has been generalized to account for all kine- 
matic regions appearing in hadronic collisions. It centers around a Lorentz invariant generaliza- 
tion of the definition of the dipole evolution variables. Splitting functions have been derived for 
the new dipole types. Together with the well-known radiation pattern of pure final-state colour 
dipoles, their utilization in a complete shower algorithm has been presented to describe soft and 
collinear multiple parton emission. The feasibility of the approach has been shown through its 
successful application to electron-positron annihilation into hadrons, inclusive Drell-Yan pair 
production and inclusive QCD jet production at hadron colliders. All comparisons deliver en- 
couraging results in good agreement with other models and with experimental data. It is worth 
to mention that the feature of generating broader pseudo-rapidity spectra - often mentioned in 
connection with colour-dipole evolution according to ARIADNE- has not been confirmed by the 
new model. Moreover, for the first time, results have been presented for the inclusive production 
of jets in hadronic collisions that have been obtained from a shower based on the colour-dipole 
approach. First evidences could be given that the model correctly accounts for colour-coherence 
effects. 

Taken together, an appealing picture of dipole cascading has been achieved. Future work will 
concern the full incorporation of gluon splittings in the initial and final state, and, the generaliza- 
tion to finite quark masses. In addition, a merging with multi-leg tree-level matrix elements for 
additional QCD radiation will be addressed, and, a matching with full NLO QCD calculations 
shall be studied. 
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